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I. INTRODUCTION 


A. BACKGROUND 

One of today’s major environmental concerns is the depletion of the 
earth’s protective ozone layer. In 1987, an international conference was 
held in Montreal, Canada to address the problems caused by 
Chlorofluorocarbons (CFCs) to the earth’s ozone layer. CFCs are manmade 
chemicals of chlorine, fluorine, and carbon and are unique in that they 
have a combination of desirable properties: low in toxicity, non- 
flammable, non-corrosive, non-explosive, extremely stable and compatible 
with many other materials. This extreme stability is what causes problems 
to the ozone layer due to the fact that CFCs only break down in the upper 
atmosphere when subjected to intense ultraviolet radiation. This break 
down produces chlorine which has been linked to the depletion of the 
earth’s ozone layer. 

In September 1987, 24 nations representing the United Nations 
Environment Program (UNEP) met and signed the Montreal Protocol. They 
discussed the substances that deplete the ozone layer [Ref.1] and called 
for a near-term freeze on the production and consumption of these 
substances. The agreement required production of these chemicals to be 
cut back to 1986 levels followed by a two-phased reduction culminating in 
cutbacks of 50% by mid-1998; this came into effect on July 1, 1989. In 
1990, a progress meeting was held in London where UNEP delegates agreed to 


completely phase out all CFCs by the year 2000 [Ref. 2]. In the spring of 


1992, President Bush pushed up the complete phase out od CFC’s by the year 
1395: 

The U.S. Navy uses a number of different CFCs (designated by "R" for 
refigerants) for various refrigeration and air conditioning (AC) needs. 
Presently, the U.S. Navy has approximately 1850 shipboard AC plants using 
both R-12 (in reciprocating compressor) and R-114 (in centrifugal 
compressor) plants. To comply with the Montreal Protocol and U.S. 
legislation, the Mechanical Systems Branch (Code 2772) at the Naval 
Surface Warfare Center (NSWC) is pursuing research in the elimination of 
shipboard use of CFCs. As mentioned by Chilman [Ref. 3] this research is 
to be completed in three phases: 

l. To identify in the short term suitable alternative ozone-safe 
chemicals to replace R-114 and R-12. To accomplish this task, the heat 
transfer characteristics must be similar to the existing refrigerants in 
place and hence the need for a database exists for current refrigerants 
(R-114 and R-12) so that they can be compared to the new proposed 
refrigerants (HFC-124 and HFC-134A respectively). 

2. In the longer term, to research, develop, and test substitute 
chemicals and alternative technologies to replace existing CFC uses. 

3. To implement new cooling system technologies into the fleet which 
do not depend on CFCs or their replacement. 

This thesis is a continuation of the previous work at NPS and 
supplements NSWC’s research on alternatives to CFCs by establishing 
baseline nucleate pool boiling data of pure R-114 and R-114/oil mixtures 
from a small bundle of enhanced tubes representing a section of a flooded 


evaporator. Emphasis is placed on the natural convection and boiling 


regions, hystersis phenomena, and analysis of various oil concentrations 


on the overall heat transfer performance. 


B. OBJECTIVES 

The objectives of this thesis are as follows: 

1. Understand in greater detail both the convection and nucleate pool 
boiling phenomena and hysteresis effects within a small Turbo-B tube 
bundle. 

2. Obtain data using a Turbo-B tube bundle for increasing and 
decreasing heat flux for R-114/oil mixture with oil concentrations of 0, 
1, 2, 3, 6, and 10 percent. 

3. Compare data with earlier studies at the Naval Postgraduate School 


using R-114/oil mixtures for other enhanced tube surfaces. 


II. LITERATURE SURVEY 


As SMOOTH TUBE BUNDLES 

In recent years significant progress has been made in understanding 
nucleate boiling heat transfer phenomena on the shell side of flooded 
evaporators. Extensive work on smooth tube bundles has been reported by 
Cornwell (Leong and Cornwell [Ref. 4], Cornwell et al. [Ref. 5], Cornwell 
and Scoones [Ref. 6], Cornwell [Ref. 7]). Cornwell and Schuller [Ref. 8] 
conducted a photographic study of boiling R-113 in a smooth tube bundle at 
one atmosphere. One of their conclusions was that bubbles leaving the 
lower tubes in the bundle impacted and caused a sliding motion around the 
upper tubes. Cornwell and Schuller observed the two-phase flow patterns 
and deduced that sliding bubbles from lower tubes on upper tubes could 
account for significant heat transfer in the top part of the bundle. 
Cornwell [Ref. 7] later found that in the nucleate boiling region, sliding 
bubbles and liquid forced convection could account for all the heat 
transfer in the top of the bundle. 

The influence of tube position within a bundle of smooth tubes using 
R-11, R-12, R-22 and R-113 has been studied extensively by Wallner [Ref. 
8], Fujita et al. [Ref. 9], Chan and Shourki [Ref. 10], Rebrov et al. 
[Ref. 11], and Marto and Anderson [Ref. 12]. Using both in-line and 
Staggered tube arrangements with various tube pitch-to-diameter ratios 
between 1.2 and 2.0, their work verified that the influence of the lower 


tubes in a bundle can significantly increase the heat transfer performance 


of upper tubes at low heat fluxes due to two-phase convective effects. At 
high heat fluxes (typically > 50 kW/m?) in the fully developed nucleate 
boiling region, the data for all the tubes merged onto a single curve. 
This is representative of a single smooth tube and shows that there is no 
‘bundle effect’ (ie. no improvement over a single tube under similar 
conditions) in the high heat flux (nucleate boiling) region. 

Chan and Shoukri [Ref. 10] studied the boiling characteristics of a 
smooth in-line tube bundle in R-113. They concluded that at lower heat 
fluxes, the heat transfer process is strongly influenced by two-phase 
convection effects, resulting in higher heat transfer coefficients on the 
upper tubes. At high heat fluxes, however, they found that the dominant 
mode of heat transfer was nucleate boiling from the upper tubes and that 
convective effects from below were insignificant. At these high fluxes, 
the bundle performance was similar to the trends of a single tube in a 
single tube apparatus. Fujita et al. [Ref. 9] also found that the heat 
transfer at low heat fluxes using a smooth tube bundle in R-113 was 
enhanced by convection induced by rising bubbles (ie. a steady increase in 
performance of the upper tubes as additional lower tubes were activated). 
They attributed this enhancement to the "positive bundle effect". A high 
heat fluxes, this enhancement disappeared. 

Anderson [Ref. 13] found similar effects as above for a smooth tube 
bundle in pure R-114. Furthermore, he reported that the presence of up to 
3% oil (by mass) actually improved the heat transfer performance. This is 
similar to data reported for a single smooth tube by Wanniarachchi et al. 
[Ref. 14]. Furthermore, at an oil concentration of 10%, only a slight 


degradation in the heat transfer (compared to pure R-114) was found. He 


obtained a maximum heat transfer performance for the bundle at an oil 


concentration of 2%. ae 


B. ENHANCED TUBE BUNDLES 

Much less work has been done on enhanced tubes (Enhanced means any 
surface that is not smooth). However, the effects at high and low heat 
fluxes mentioned above are similar to those obtained for finned tube 
bundles by Yilmaz and Palen [Ref. 15], Muller [Ref. 16], and Hahne and 
Muller [Ref. 17]. Stephen and Mitrovic [Ref. 18] looked at R-12 and R-114 
boiling from a GEWA-T tube bundle. Apart from the magnitude of the heat 
transfer coefficient varying with fluid, the trends are very similar to 
those mentioned above for smooth and finned tube bundles. 

For porous coated surfaces, Czikk et al. [Ref. 19] found no ‘bundle 
effect’ over a wide range of heat flux (1-100 kW/m*) and the bundle data 
agreed closely with single tube data. Arai et al. [Ref. 20] found that 
the ‘bundle effect’ for a Thermoexcel tube bundle was smaller than that 
found for a smooth or finned tube bundle. However as before, any ‘bundle 
effect’ was eliminated at high heat fluxes where the data for all the 
tubes agreed closely with single Thermoexcel tube results. These effects 
are similar to those found by Czikk et al. [Ref. 19] for the porous 
coated. 

Chilman [Ref. 3] reported experiments with a Turbo-B tube bundle in 
pure R-113, conducting both increasing and decreasing heat flux tests. He 
concluded that in the natural convection region, heated lower tubes do not 


have any influence on the heat transfer from upper tubes. Also, Chilman 


reported the presence of heated lower tubes within a bundle reduced the 
incipient boiling point. 

Stephan and Mitrovic [Ref. 18] reported the influence of oil on the 
boiling heat-transfer coefficient of R-12 using a T-shaped finned tube 
(Gewa-T) bundle. They reported that the ratio of oil to no oil heat- 
transfer coefficients decreased with the mass fraction of oil for all 
except the highest heat flux (22 kW/m’) where an increase in heat transfer 
was noted for oil concentrations between 1 and 6% . They concluded that 
the influence of oil on heat transfer was mainly due to the thermal 
properties of the specific oil used in the experiments and its interaction 
with the refrigerant. 

Schlager et al. [{Ref. 21] summarized the influence of oil on 
refrigerant in pool boiling. They stated that under certain conditions 
(typically low pressure and high heat flux), the heat transfer coefficient 
increased at low oil concentration. Stephan [Ref. 22] first pointed this 
out and attributed the phenomenon to foaming. Burkhardt and Hahne [Ref. 
23] for a finned tube bundle found that the maximum heat transfer 
coefficient, which was 10% to 15% above the oil-free value, occurred at a 
concentration of about 4%. 

Heimbach [Ref. 24] conducted experiments with R-12/oil mixtures on a 
finned tube bundle. He reported that the presence of up to 2% oil, did 
not affect the heat transfer performance significantly. However at higher 
concentrations (3% to 7%), an increase in the heat transfer was observed. 
He also attributed this to foaming and postulated that changes in the 


properties of the mixture might facilitate the formation of bubbles. 


Anderson [Ref. 13] and Akcasayar [Ref. 25] conducted experiments with 
finned (19 fpi) and High Flux (porous coated) tube bundles in pure R-114 
and R-114/oil mixtures in the same apparatus. Anderson reported a maximum 
heat transfer performance at an oil concentration of 3% for the finned 
tube bundle. Akcasayar also reported that the finned tube bundle 
performance increased 1.65 times with 3% oil concentration (compared with 
pure R-114) at the maximum heat flux level. For 6% and 10% oil 
concentrations, the performance of the bundle, when compared with lower 
oil concentrations, decreased. This was especially significant at a 10% 
oil concentration. Compared with the finned tube bundle, the High Flux 
tube bundle had a 1.5 times better heat transfer performance at a heat 
flux of 30 kW/m* in pure R-114. However, these performance ratios 
decreased with increased oil concentrations such that the finned bundle 
outperformed the High Flux bundle at 6% oil concentration. This was 
especially true at the highest heat fluxes where the High Flux bundle 


performance was not much better than a smooth tube bundle of similar size. 


III. EXPERIMENTAL APPARATUS 


A. TEST APPARATUS OVERVIEW 

The experimental apparatus including the auxiliary equipment and the 
evaporator/condenser is shown in Figure l. The following is only a 
general description of the whole experimental apparatus. A more detailed 
look at the condenser and evaporator is provided in section C. Further 
information about the apparatus is provided by Murphy [Ref. 26] and 
Anderson [Ref. 13}. 

The experimental apparatus is essentially made up of three closed 
loops. The first loop consists of an 8 ton refrigeration unit located 
outside the laboratory which is used to cool an ethylene glycol/water 
mixture. The second loop is the ethylene glycol/water mixture flowing 
through the condenser. This mixture is contained within a large sump 
within the laboratory. The flow rate through the condenser is delivered 
by two pumps which can be operated independently or together; this coolant 
mixture condenses the refrigerant vapor in the condenser and maintains 
system pressure and temperature. Pump number one provides coolant flow 
through the four test condenser tubes as well as to one of the auxiliary 
condenser coils (bottom coil). Pump number two provides coolant through 
a manifold which distributes the coolant to the remaining four auxiliary 
condenser coils within the condenser. The third loop is the evaporator 


and condenser itself designed for reflux operation. The vapor generated 


in the evaporator flows upward and condenses in the condenser; the 


condensate then returns to the evaporator via gravity. 


Ba AUXILIARY EQUIPMENT 
i. 28 kW Refrigeration Unit 
This unit was used to cool a 1.8 mm reservoir sump of an 
ethylene glycol/water mixture (coolant) to the desired temperature needed 
to condense the refrigerant vapor. For R-114, the temperature control was 
set to maintain the sump at -15 °C. The refrigeration unit had a cooling 
capacity of 28 kW (8 tons). 
2: Ethylene Glycol/Water Mixture 
The coolant used was a 54% (by weight) ethylene glycol/water 
mixture. This refrigerated mixture was used to control the system 
pressure and temperature by circulating coolant through the auxiliary 
condenser coils and/or condenser test tubes at different flow rates. 
oe Pumps 
Two pumps were available to circulate the coolant from the sump 
through the condenser. Pump number one fed four test condenser tubes and 
one of the auxiliary condenser coils. Pump number two fed the other four 
auxiliary condenser coils. Pump number one was the primary pump used at 
low heat fluxes; pump number two was started (as necessary) at high heat 
fluxes to maintain the desired saturation pressure in the evaporator. 
4, Flowmeters 
Five calibrated float-type flowmeters, connected to pump number 
one, were used to measure the flow rates passing through the four test 


condenser tubes and one auxiliary condenser coil. One additional 
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flowmeter (connected to pump number two) was used to measure the total 
flow to the remaining four auxiliary condenser coils. Each of the five 
auxiliary condenser coils had a globe valve to regulate (or shut off) 
flow as desired. For the four test condenser tubes, the coolant flow was 


regulated by a flowmeter valve. 


C. EVAPORATOR/CONDENSER 

An overall view of the evaporator and condenser is shown in Figure 2. 
The evaporator was designed to simulate a small portion of a refrigerant 
flooded evaporator. Front and side views of the evaporator are shown in 
Figures 3 and 4. It was fabricated from stainless steel plate and formed 
into a short cylinder, 610 mm in diameter and 241 mm long. Electrically- 
heated tubes were cantilever—mounted from the back wall of the evaporator 
to permit viewing along the axis of the tubes through the lower of two 
viewing windows. A plexiglas plate was attached to the front of the tube 
bundle to ensure tube alignment during experiments. Each viewing window 
had a layer of glass and plexiglas, the glass being used on the 
refrigerant side in order to prevent surface cracking of the plexiglas. 
The plexiglas gave the glass added strength and served as a safety barrier 
in case the glass broke during pressurization. 

The electric power can be applied separately to each set of heaters 
using a STACO 240 V, 23.5 kVA rheostat controller shown in Figure 5. 
Also, the desired number of instrumented tubes, active tubes, simulation 
or auxiliary heaters can be individually activated by using circuit 
breakers. The five simulation heaters, each with a maximum rating of 4 


kW, were mounted below the test bundle in order to simulate 15 additional 


ll 


tube rows in a larger bundle and to provide an inlet vapor quality into 
the bottom of the test bundle as suggested by Webb [Ref. 27]. The liquid 
pool was maintained at 2.2 °C (corresponding to a saturation pressure of 
l atmosphere) by passing coolant through the condenser. 

Figure 6 is a schematic sectional view of the evaporator that shows 
the four kinds of heated tubes installed in the evaporator. These were 
instrumented, active, auxiliary, and simulation. For this study only the 
instrumented, active, and simulation heaters were used; the auxiliary 
heaters are needed for experiments either at higher pressures or for other 
refrigerants which have a higher normal boiling point (such as R-113). 
Table 1 gives the power rating for these heaters and the number used in 
the evaporator. 

The tube bundle itself consists of instrumented, active, and dummy 
tubes. The location of each tube is represented by the respective letter 
I, A, and D as shown in Figure 6. The test bundle consists of two types of 
heated tubes: 12 active tubes (marked "A") which contained 1 kW cartridge 
heaters, and 5 instrumented tubes (marked "I") which, in addition to the 
1 kW cartridge heaters, contained six wall thermocouples each. 

In measuring boiling heat transfer coefficients, great care must be 
exercised with the cartridge heater and temperature measuring 
instrumentation to ensure good accuracy. The instrumented test tubes were 
fabricated using the same method as that used by Hahne and Muller (Ref. 
17] and Wanniarachchi et al. [Ref. 14]. The exact procedure can be found 
in Eraydin (Ref. 28]. Figure / is a cross-sectional sketch of an 
instrumented tube, showing the construction details and the location of 


the wall thermocouples. The thermocouples were embedded in the wall at 
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different circumferential and longitudinal positions along the heated 
section of the tube shown in Figure /. 

The five instrumented tubes were located along the centerline of the 
tube bundle, forming a vertical in-line column. All the instrumented and 
active tubes were Turbo-B tubes made by Wolverine Tube Co. (see Section 
E). These tubes were nominally 15.9 mm in outside diameter and were 
arranged in an equilateral triangular pitch (ie. centerline—to-—centerline 
spacing) of 19.1 mm, giving a pitch-to-diameter ratio of 1.35. The 
thickness of the Turbo-B enhancement was 0.85 mm giving a diameter to the 
base of the enhancement of 14.2 mn. 

The bundle also contained a number of unheated dummy smooth tubes 
(marked "D") that were used to guide the two phase mixture through the 
bundle. The dummy tubes were made from commercially available 15.9 mm OD 
smooth copper tubing. Two vertical baffle plates made of aluminum were 
used on either side of the bundle to restrict circulation into and out of 
the bundle at the sides. A dummy rack (Figure 8) consisting of 12 solid 
rods made of aluminum (15.9 mm OD and spaced 19.1 mm from centerline-to-— 
centerline) was placed below the tube bundle. This rack had a triangular 
pitch arrangement with vapor retaining plates on the sides and was 
designed for two purposes: to collect all rising two phase flow generated 
by the simulation heaters and direct it into the test bundle and to 
simulate vapor passing through a large bundle before reaching the 
instrumented tubes. A small open space (approximately 5 mm in height) was 
left between the bundle and dummy tube rack. This space allowed some 
refrigerant from below to enter the bundle and replace the vapor being 


generated in the bundle. However, there was also a space below the dummy 
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rack that allowed the majority of the circulation to occur. Thus, 
liquid/vapor circulation was vertically upward over the five instrumented 
test tubes with no net horizontal component. Most of the liquid—vapor 
mixture after passing through the bundle was separated when it reached the 
pool surface. However, due to the strong circulation patterns set up 
within the liquid pool, some vapor bubbles remained trapped in the liquid 
and circulated around the pool. 

The condenser included four test tubes (each 1.219 m in length and 
15.9 mm OD) in a vertical in-line column and five auxiliary copper coils. 
For the boiling experiments, these tubes were used to regulate the 
pressure and temperature in the evaporator. The condenser was designed to 
permit independent condensation studies of small in-line tube bundles, 
using the evaporator as a source of vapor. Details of the condenser can 


be found in Mazzone [Ref. 29]. 


Di DATA ACQUISITION SYSTEM/INSTRUMENTATION 

As described by Akcasayar [Ref. 25], a Hewlett Packard HP—3497A Data 
Acquisition System, HP-9125 computer and HP—/01 printer were used for data 
acquisition, data reduction and data printing respectively. Although an 
HP-9826 computer and HP-/7470A plotter can be used for final graph 
printing, a Macintosh computer (using Criketgraph) was utilized. As 
described by Anderson [Ref. 13], type-T copper-constantan thermocouples 
were used to make temperature measurements on the HP-3497A using a relay 
multiplexer assembly equipped with thermocouple compensation. A 20— 
channel relay multiplexer card was used to measure the voltage output from 


both the voltage and current sensors. The voltage measurements were taken 


14 


from separate sensors that measured the voltage going to the tube bundle, 
the simulation heaters and the auxiliary heaters. The total current 
going through the auxiliary and simulation heaters was measured using an 
American Aerospace Control (ACC) current sensor. The current to each 
instrumented tube heater was measured using five identical current 
sensors. The voltage supplied to the other active tubes was also 
measured, but the current of each active tube was not. Instead, the total 
current for a pair of active tubes was measured, and this was sufficient 
since these tubes each had the same power output (1000 W) as the 
instrumented test tube heaters and there was no apparent reason to monitor 
each active tube individually. Computer channel assignments for data 


acquisition and array assignments are given in Table 2. 


EE. GEOMETRY OF TURBO-B TUBE 

The Turbo—-B tube, manufactured by Wolverine Tube Inc., contained an 
enhanced surface geometry. The exterior boiling enhancement is made by 
raising low integral fins, cutting diagonally across these fins, and then 
rolling the fins to compress them to form mushroom—like pedestals [Ref. 
30]. This process forms numerous re-entrant passageways. Figure 9 shows 
the surface of the tube at 25 times its actual size. The tube is 


currently available in copper, cupro-nickel, and low carbon steel. 


The relative dimensions of the tube used in this study are as follows: 
Tube material — Copper 
Nominal Outside diameter = 15.9 mm 


Enhanced surface length = 203.2 mm 
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Thickness of Enhancement = 0.85 mm 


Diameter to Base of Enhancement = 14.2 mm 
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Figure 1. Schematic View of the Experimental Apparatus 
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Figure 2. Evaporator/Condenser Schematic 
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Figure 3. Front View of Evaporator 
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Figure 4. Side View of Evaporator 
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Figure 5. Photograph of 208 V, 75 A, Variable Transformers 
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Sectional View of Evaporator Showing Tube Bundle 


Figure 6. 
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Thermocouple Locations on an Instrumented Boiling Tube and 


Tube Section View 
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Photograph of Dummy Rack 


Figure 8. 





Table 1. EVAPORATOR HEATERS 


Heater Type 


instrumented Tube Heaters | 5 | 
- 43 










Power Rating per Heater 







Active Tube Heaters 
Auxiliary Heaters 


Simulation Heaters 





Table 2. COMPUTER/DATA ACQUISITION 
ASSIGNMENT 


= 


Amperage Sensor Description 


mr a A 










tube? Sd SO *dYSSCAm (0) 
[ube 28d Amt 
[Tube 3 SS—~sCSSS—tSdSSC Amt) 
[Tube 4 SSS Sid Am) 
[ube 5 S—S—sSSCS~—i Sd Am) 


|__Active Heater Group 1 | 35 | Amp (5) 
|__ Active Heater Group 2. | 36 | Amp (6) 
|___Active Heater Group 4 | 38 | Amp (8) 
|__ Active Heater Group 5 | 39] mp9) 
__ Auxiliary Heaters | 25] Amp(10) 
|___Simulation Heaters | 26] Amp (14) 








|____ Simulation Heaters | 28] Voit(t) 
|____ Auxiliary Heaters | 29] Voit (2) 
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Table 2. COMPUTER/DATA ACQUISITION ASSIGNMENT (CONT. ) (cont. ) 


| Vapor 1-Top of Condenser | 00s | (0) 

Vapor 3-Top of Evaporator | 02 #| (2) | 
| _Liquid 1-Top of bundle =| 03 | (3) 
| _Liquid 2-Top of bundle | 04 | (4) 
“Tiquid S-Bottom of bundle] 05 | 15) 
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TT ee ee a ae 
SSS EO ae a a 
bem Noms 4 9) 
po Tube No.5 TNO) 
eel Nore 4s 1) 
ee re ae a aa ce 
ber Nocmmn ten 480 fT (14) 
era Nora 4 5) 
ES a ee ae 
eA Ge a al a cr ae 
TE a a 
| Tube Now} ST (19) 
| Tube 3,No.3 0 T5420) 
aT ibels Nos SS Tet) 
bora Novo nese ez) 
abe 4iNow sa ie Tes) 
ices Nore | se 2s) 
| Tube 4,No.300 TO (26) 
ers Nonsense Pez Ns) (28) 
| Tube 4.No.6 TS T(29) 
_ EES nn a a ac 
ES De eee ee ee a 
oer Norse | es P32) 


Tube 5,No. 5 ae T(34) 
Tube 5,No. 6 pS (35) 
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IV. EXPERIMENTAL PROCEDURES 


A. REMOVAL OF THE TUBE BUNDLE AND BUNDLE DISASSEMBLY 

Before starting the removal of the tube bundle from the evaporator, 
the front glass viewing windows were carefully removed. Next, all 
thermocouple wires and tube heater electrical connections’ were 
disconnected. After this was completed, the nuts securing the backing 
plate and support block were removed, and the tube bundle was taken out 
from the back of the evaporator. 

When the bundle needed to be disassembled, it was ensured that there 
was a clean working surface. The first task was to remove the plexiglas 
plate attached to one end of the aluminum baffle plates (ie. the front of 
the bundle assembly) by four screws. The ten screws on the side of each 
aluminum baffle plate were then removed (these were attached to the dummy 
tubes down each side of the bundle). The aluminum plates were then pulled 
off the bundle. The four corner dummy tubes (two top and two bottom) 
remained attached to the tube bundle support block as they were 
countersunk into the block. The six outer smooth tubes (three per side) 
could be easily pulled from the aluminum plates as they were attached only 
by the screws already removed. The other ten smooth tubes were then 
unscrewed from the tube bundle support block as seen in Figure 10. The 
smooth tubes were engraved to ensure proper identification during 
reassembly. With these tubes and aluminum baffle plates removed, only the 


instrumented and active enhanced heater tubes remained. These tubes were 
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removed from the support block by loosening the outer O-ring compression 
plate, disconnecting the active heater tube wired pairs, and pulling the 
tubes from the block. Reassembly of the tube bundle is done by reversal 


of this procedure. 


Be SYSTEM CLEAN-UP 

If the system had been previously operated with refrigeration oil (or 
contaminated from some other source) it had to be thoroughly cleaned. To 
accomplish this, the entire apparatus had to be taken apart and cleaned in 
the following manner. 

After removal of the refrigerant (R-113 by directly draining into 5 
gallon drums via drain valve R-5 (see Figure 2 in Chapter 3) at the bottom 
of the evaporator or R-114 by boiling off into the storage tank by opening 
R-1l and R-/) and with the system at atmospheric pressure, all electrical 
connections to the bundle were disconnected and the front viewing glass 
windows were removed. An electric fan was used for safety to ensure 
proper ventilation. The tube bundle was then removed and disassembled as 
described in section A; the dummy tube rack was also removed. 

Having removed the tubes from the tube bundle, they were individually 
washed with warm water, rinsed and then wiped down with acetone. The 
smooth tubes were cleaned with Copper Brite (a commercial copper cleaning 
product) to remove any oxidation. They were also wiped down with warm 
water and then with acetone. The same procedure was followed for the 
Turbo-B tubes except they were not cleaned with Copper Brite for fear of 
clogging the channels. During the cleaning process, a soft bristled 


toothbrush was used to ensure the enhanced surface was cleaned properly, 
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exercising care not to interfere with the tube surface. The evaporator 


shell was cleaned in a similar manner, using warm water and acetone. 


GC. INSTALLATION OF THE TUBE BUNDLE 

Once the tube bundle had been cleaned and reassembled (see section A), 
and before tightening the backing plate nuts, the whole assembly was 
carefully guided back into the evaporator section, ensuring the plexiglas 
viewing cover of the tube bundle was not damaged. After the bundle was in 
position, it was ensured that the dummy tube rack was properly positioned 
below the bundle and that the vapor thermocouple positions were still 1.75 
cm above the bundle. Then, all the nuts were tightened equally on opposite 
sides to give equal compression on the gasket. To replace the front 
window, very small, equal torques (using a torque wrench) were applied 
circumferentially to each nut on the outer ring support in turn. After 
the window was in place, each tube (which extended through the outer O- 
ring compression plate) was lightly tapped forward so as to touch the 
front-—viewing window. The backing plate was then tightened and the 
individual tube O-rings compressed, providing a good seal for the system. 
The compression plate had grooves for the tube O-rings to sit in to help 


with proper alignment and ensure a good seal. 


D. SYSTEM LEAKAGE TEST 

After the system was isolated from the atmosphere and system integrity 
was restored, a Seargent Welch 10 SCFM vacuum pump was connected to the 
apparatus (via valves R-l and R-8) and the pressure taken down to 25 inHg 


vacuum. Valves R-1l and R-8 were then secured and the system was left 
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untouched for at least 24 hours to see if there was any air leakage in. 
If there was significant leakage (>l inHg over 24 hours), then the vacuum 
was broken by cracking open valve R-2 slowly (this ensured that no 
moisture entered the system). The system was pressurized (with air) to 15 
psig through valve R-2. Large leaks could then be detected by simply 
listening to the air issuing from the system; small leaks were detected by 
spraying a soapy water solution to all surfaces where leaks were most 
likely to occur (front viewing glass gaskets, backing plate gasket, all 
fittings/valves coming off the condenser/evaporator, O-ring seals of the 
bundle tubes etc). Extreme care must be taken to ensure no moisture 
enters ene inside of the heated tubes where the heater wires protrude. 
After all leaks were detected and corrected, the system was again 
subjected to a vacuum for another 24 hour period. If the vacuum held, 
then the system was ready to receive refrigerant. If not, the above leak 


correction test was repeated. 


E. REFRIGERANT 
i Fill 
a. From System Storage Tank 

A refrigerant storage tank was used to store R-114 during 
modification/repairs to the system. The storage tank prevented discharge 
of the R-114 into the atmosphere and made the experimentation less costly. 
To fill the evaporator with R-114 from the storage tank, the ethylene 
glycol/water coolant temperature was first reduced to -15 °C. The system 
pressure was then maintained below the storage tank pressure (vapor 


pressure of R-114 at 20 °C is approximately 15 psig) by circulation of the 
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coolant through the condenser test tubes and auxiliary coils. Valves R-6 
and R-4 were then opened to draw the R-114 from the storage tank to the 
evaporator. The amount of refrigerant that was transferred was controlled 
by throttling valve R-6 to obtain the desired level in the evaporator. If 
required, additional R-114 could be transferred from a 68 kg storage 
cylinder to the system using valve R-2 (see section 1.b). 
b. From Refrigerant Storage Cylinder 
To fill the apparatus from the 68 kg storage cylinder, the 
system pressure was reduced in the same way as above. A hose assembly 
containing a Drierite gas purifier was connected between the storage 
cylinder and valve R-2. A_ gas purifier was used not only to remove all 
impurities, but also to remove any water from the refrigerant. Once in 
place, both the storage cylinder valve and R-2 were opened until the 
desired refrigerant level was reached in the evaporator. 
7 is Removal to the Storage Tank 
For tube replacement, system maintenance or system clean up, the 
R-114 was transferred to the storage tank. The ethylene glycol/water 
coolant temperature flowing through the storage tank was cooled to -15 °%; 
valves R-/ and R-8 were opened and the vacuum pump was turned on to put 
the storage tank under vacuum. Once the storage tank was under a 20 inHg 
vacuum, valve R-8 was shut and valve R-1 was opened. Next R-114 was slowly 
boiled off to the storage tank by using the tube bundle, simulation and 
auxiliary heaters at a heat flux of 600 kW/m* (slow boiling is important 
to ensure minimum transfer of oil). As the refrigerant level decreased, 


individual heaters were turned off to ensure none were uncovered. Once 


a2 


below the level of the heaters, the final few cm of R-114 was boiled off 
using heat from the atmosphere. Once all of the R-114 was transferred, 


valves R-l and R-/ were shut. 


Be. OPERATION 
1a System Startup, Securing and Emergency Procedures 
See Appendix D 
2. Normal Operation 

The evaporator was filled with R-114 to a level of approximately 
10 cm above the top tubes in the bundle. Prior to operating the system, 
the 8 ton refrigeration unit was run for approximately an hour to reduce 
the ethylene glycol/water coolant in the sump to a temperature of -15 °C. 
The pressure in the evaporator/condenser was usually 12 to 15 psig if the 
system had been secured overnight. As the sump was brought to 
temperature, the data acquisition system and computer were turned on. 
This allowed the temperature in the system to be monitored during cool- 
down to saturation conditions. With this and pump number one running, 
one auxiliary condenser coil and the four condenser test tubes were used 
to bring the pool down to a subcooled condition (for R-114, approximately 
1 °C on all _ three pool thermocouples). Subcooling of the refrigerant was 
done to ensure the pool had an evenly distributed temperature prior to 
Starting a run. After reaching this subcooled condition, all coolant 
supply to the condenser was secured. The pool was then allowed to ‘heat 
up’ by conduction from the surroundings. Once a saturation temperature of 
2.2 °C was reached, the instrumented tube(s) (and simulation heaters for 


test 7) was/were switched on and set to the desired heat flux value. 
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This lengthy procedure was done to prevent the tubes from prematurely 
nucleating. The heat flux of the instrumented tubes was then slowly 
increased at desired intervals by adjusting the rheostat. For increasing 
heat flux, the data was taken with very small heat flux increments (every 
1000 kW/m?) , waiting at least 5 minutes to attain steady state conditions 
at each heat flux. At all regions of the boiling curve (and especially 
near the onset of nucleate boiling), two readings were taken at each heat 
flux to ensure accuracy. The bundle was continuously visually monitored 
through the observation windows. Figure 11 shows the tube bundle 
arrangements used during the experimentation. Test one was with only one 
instrumented tube turned on at any position within the bundle. Test two 
was with instrumented tubes one and two active in the bundle. Test three 
was with instrumented tubes one, two and three active in the bundle. Test 
four was with instrumented tubes one, two, three, and four active in the 
bundle. Test five was with all five instrumented tubes active in the 
bundle. Test six was all five instrumented tubes plus all five pairs of 
active enhanced heater tubes active in the bundle. Test seven was the 
same aS test six with the addition of all five simulation heaters active. 
For each data set, the five simulation heaters had the same heat flux as 


the tubes within the bundle. 


G. OIL ADDITION 

During the bundle experiments, successive amounts of York-C oil were 
added into the evaporator. Since the weight of the refrigerant in the 
evaporator was 60.3 kg, the amount of oil corresponding to 1% by weight 


was measured as 6/0 ml, 2% 1340 ml etc. The oil was syphoned into the 
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evaporator via a funnel/hose connection through valve R-3 by reducing the 
pressure in the evaporator to less than 15 inHg vacuum. Ensuring that no 
air entered the system, valve R-3 was promptly shut when the desired 


amount of oil had been added. 


H. DATA REDUCTION PROCEDURES 

The data reduction program "DRP4RH" was used during the experiments 
for processing the data collected ( see Appendix E for listing). The 
program was written in HP Basic 3.01 and run on an HP-9000 series 
computer. The characteristics and capabilities of this software are 
similar to those provided by Anderson ([Ref. 13]. The following 
modifications were made: 
1. Correction for pool height by Chilman [Ref. 3] 
2. Installation of new thermocouple at the bottom of the liquid pool 
(bundle inlet temperature) by Chilman [Ref. 3] 
3. The ability to obtain data from one instrumented tube at any position 
in the bundle. 
4. The natural convection correlation of Churchill and Chu [Ref. 31] for 
a single horizontal cylinder in an ‘infinite’ liquid pool was added for 


comparison with experimental data. 
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Figure 10. Photograph of Tube Bundle support Block 
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Figure 11. Tube Bundle Arrangements used During Experimentation 
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V. RESULTS AND DISCUSSION 


Ay INTRODUCTION 

The results are presented in four sections with subs-sections as 
appropriate. The first section discusses the preliminary experiments 
which led to modifying the experimental start-up procedure to include 
subcooling. The second section discusses the natural convection effects, 
nucleate pool boiling phenomena and hysteresis effects within the tube 
bundle in pure R-114. The third section discusses similar phenomena for 
R-114/oil mixtures and their effects on the above. The fourth section 
shows comparisons of data taken during this thesis with previously 
obtained data at the Naval Postgraduate School. 

A list of data files taken during this investigation may be found in 
Appendix A. All data files used in this thesis use the following filename 
sequence. Each file is composed of five sets of alpha-numeric characters 
used to describe the experiment. 

First set (2 char.) Tube Type TB (Turbo-B) 


Second set (1 char.) Heat Flux I (Increasing) 
D (Decreasing) 


Third set (2 char.) Oil Percent 00 (0%) 
O01 (1%) 
02 (2%) 
03 (3%) 
06 (6%) 
10 (10%) 


Fourth set (2 char.) Test Type Ol (test 1) 
O02 (test 2) 
03 (test 3) 
04 (test 4) 
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O5 (test 5) 
O06 (test 6) 
O7 (test 7) 

Fifth Set (1 char.) Additional tests A-Z (If conducted) 

To give an example, the filename TBIO107 means "Turbo-B tube bundle, 
increasing heat flux, 1% R-114/oil mixture and test number 7". If more 
detail is desired about a specific data set, see Appendix A. All plot 
filenames are similar to data filenames except they start with the letter 
"Pp". The test numbers are shown in Figure 11 in Chapter IV. 

All graphs are plotted showing heat flux (W/m*) along the ordinate (y 
axis) as a function of wall superheat (K) along the abscissa (x axis). 
The wall superheat is defined as the difference between the ‘corrected’ 
average tube wall temperature (ie. having accounted for depth of 
thermocouple burial) and the local liquid saturation temperature 
(corrected for hydrostatic head within the bundle). The heat flux was 
corrected to account for the heat lost through the unheated tube ends. 
The heat flux was varied from 600 to 100,000 W/m* for increasing heat flux. 
To ensure greater detail at the point of incipience, the heat flux was 
increased in small steps; these settings varied from test to test. The 
heat flux values for decreasing experiments were taken at prescribed 
settings for easy comparison with past experiments and future reference. 
Diese heat flux settings were 1xl10°, 7.5x10*, 5x10°, 3x10*, 2x10*, 1.5x10*, 
1x10*, 7x10°,4x10°, 2x10°, 1x10° W/m?. Approximately 30-40 data points were 
taken for each increasing heat flux run and 20-25 points were taken for 


each decreasing heat flux run. 
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B. PRELIMINARY EXPERIMENTS 

After the pure R-113 was removed following Chilman’s [Ref. 3] 
experiments and the apparatus and system cleaned, pure R-114 was added 
from the storage cylinders. Five tests were conducted using this R-114 
and these are shown in Figure 12. The first test (TBIOOO1A) was test one 
with the top tube activated. The procedure described by Eraydin [Ref. 28] 
was followed, but the plot of test TBIOOOI1A and Eraydin’s data for test 
one (also shown in Figure 12) produced significantly different results in 
the natural convection (NC) region. The data of Eraydin show a greater 
heat transfer coefficient (lower wall superheat) than test TBIOOO1A which 
show results closer to the Churchill/Chu [Ref. 31] correlation (C/C) for 
natural convection. The only difference in the apparatus between test 
TBIOOO1LA and Eyradin'’s experiments was the addition of a third 
thermocouple at the bottom of the pool. For test TBIOOOI1A, all three pool 
temperature readings were within +/- 0.1 °C prior to recording data. For 
the data of Eraydin, only the temperature at the top of the pool could be 
checked. 

Test TBIOOO1B was thought to be a repeat of test TBIOOOIA, but for the 
bottom tube in the bundle (tube 5 only). However, upon observation and 
investigation, it was found that two tubes (tube one and tube five) were 
activated due to the way the program DRP4 was set up. Hence the data 
presented is for tube one (the top tube) with tube five (the bottom tube) 
activated as well. The program was then modified to obtain data for a 
single tube (test one) at any position within the bundle. 

Test TBIOOO1C was conducted using the bottom tube (tube five) as a 


single tube. Again following Eraydin’s [Ref. 28] procedure, partial 
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nucleation was observed immediately, explaining why the data lie well to 
the left of the Churchill/Chu correlation (ie. a higher heat transfer 
coefficient). It can be seen that test TBIOOO1C is in good agreement with 
the single tube of Eraydin (top tube). The reason for this is probably 
due to the start-up conditions, which were not carefully monitored for any 
of these tests. 

Test TBIOOO1D was conducted using only tube three and also followed 
the procedure of Eraydin. This displayed similar behavior as test 
TBIOOO1A (tube one) except that nucleation was delayed, ie. occurred at a 
higher wall superheat. 

All experiments thus far were conducted following Eraydin’s procedure 
for pure R-114. It was next decided to vary his procedure slightly. For 
test TBIOOO1E, the pool was first subcooled slightly to 1°C. This ensured 
that any nucleation sites were deactivated and that the whole pool was at 
an even temperature. The pool was then brought up to the required 
saturation temperature. However, partial nucleation was still observed 
during the run. The plot of TBIOOO1LE is similar to test TBIOOOIC and 
Eraydin’s test one. 

With all of these confusing results, it was decided to empty, clean 
and recharge the evaporator with fresh R-114. Upon boiling off the ‘old’ 
R-114, a small quantity of oil contamination was found in the bottom of 
the evaporator. This could have been either vacuum pump oil (which may 
have leaked into the apparatus) or refrigerant oil that entered the system 
with the R-114 charge (which should be minimal). A third possibility 
could be oil from previous refrigerant/oil mixture experiments. However, 


this is not likely since the system had been completely stripped and 
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cleaned twice since the last mixture tests had been conducted. Samples of 
each type of oil (vacuum pump oil, miscible refrigerant oil (York-C) and 
a sample of the contamination) were sent away for analysis and the results 
were still forthcoming at the time of writing this thesis. From the color 
of the contamination, it would appear to have originated from the vacuum 
pump. This would seem to be more logical since Chilman had experienced 
problems with the vacuum pump. 

The system, including the bundle, were cleaned thoroughly and fresh 
R-114 was added. However, a gas purifier was utilized to ensure that only 
pure, clean R-114 (with no moisture) was added (see Chapter IV for further 
details). Test TBIOOO1F then was conducted following Eraydin’s [Ref. 28] 
procedure with no subcooling and ‘ignoring’ the bottom pool thermocouple 
value. The data showed good agreement with Eraydin in the natural 
convection region, but a lower heat-transfer coefficient in the boiling 
region. A possible explanation for this lower heat transfer is the fact 
that each individual tube was fully cleaned prior to adding the new R-114 
and the surface characteristics may have been modified in some way. 

Test TBIOOO1G was conducted using the same procedure except the pool 
was initially subcooled down to 1°C for 30 minutes. This ensured the pool 
had an even temperature distribution throughout. The data were now much 
closer to the Churchill/Chu [{Ref. 31] correlation. However, some 
premature nucleation was still observed. Test TBIOOO1J repeated the above 
test with the pool temperature subcooled to 1°C for at least one hour to 
further deactivate any remaining nucleation sites within the bundle prior 
to starting experiments. The data then agreed with the Churchill/Chu 


correlation as seen in Figure 12. It became apparent that premature 
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nucleation could affect the natural convection data significantly. 
Therefore, for all subsequent tests, this same procedure was adopted with 
strict observation of the bundle to ensure no premature nucleation 


occurred. 


C. PURE R-114 TURBO-B TUBE BUNDLE EXPERIMENTS 

i Test One for Different Tube Positions 
The first set of experiments conducted were performance tests in 
pure R-114. Figure 13 shows increasing heat flux for a single tube within 
the bundle at different positions (positions 1, 3, and 5) while Figure 14 
shows corresponding data (including typical uncertainties) for decreasing 
heat flux. All three tube positions agree closely with the Churchill/Chu 
[Ref. 31] correlation in the natural convection region. The difference 
between position 1 and position 3 and 5 may be that position 1 is affected 
by the fact that the flow is ‘free’ to expand after leaving the bundle. 
This difference may also be due to wall temperature uncertainty due to 
differences in the fabrication process (see uncertainty analysis Appendix 
C); however in the natural convection region, this uncertainty is low due 
to relatively high values of wall superheat. Figure 13 also shows that 
tube position within the bundle may influence the point of incipience. 
Bergles and Rohsenow [{Ref. 32] have studied the incipient point in more 
detail. They concluded that nucleation was controlled by a nucleation 


parameter, N, given by 


(7) (Ts) (Vg) 


(py) (h¢g) 
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which creates an incipient boiling superheat given by 

(T, - Ts) = 2(N)/r. 

In the above expression, r is the local bubble radius. Calculation of the 
nucleation parameter for pure R-114 showed that as saturation pressure 
increases, the nucleation parameter decreases. Assuming that the radius 
of curvature of the forming bubbles is constant (which is reasonable for 
a Turbo-B surface which has large, regularly spaced cavities), then (T, — 
T,) also decreases and nucleation may be expected to occur earlier. This 
certainly seems to be verified in Figure 13 where the incipient point 
occurs earlier (lower wall superheat) for a lower tube ie. where there is 
an increase in the local saturation pressure. You et al. [Ref. 33] also 
showed a decrease (approx. 30%) in the average incipient superheat as 
pressure was increased from 1 to 1.5 bar for pool boiling of FC-72 on a 
single tube, offering some other experimental verification for this 
conclusion. 

Chilman [Ref. 3] conducted test one for R-113 using the top tube 
only and varied the local saturation pressure by varying the pool height 
in the evaporator. He found that the point of incipience was delayed when 
increasing the liquid pool height (ie. the hydrostatic pressure head). 
Boundary layer effects due to different liquid circulation patterns may 
have caused this delay in nucleation and more research is certainly needed 
to fully understand the influence of pressure on the point of incipience. 
One experiment that could be conducted would be to vary the pool height, 
but keep the local pressure at each tube constant by simultaneously 


varying the vapor pressure above the pool. 
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Once nucleation occurs, Figure 13 shows that the single tube 
experiments merge onto a single boiling curve. In this region, there 
appears to be no effect of hydrostatic pressure head. Figure 14 shows the 
corresponding decreasing heat flux data for a single tube at the same 
three positions within the bundle. It shows no significant influence of 
tube position in the bundle for decreasing heat flux. Note that at low 
heat flux, the experimental uncertainty in heat flux and wall superheat is 
the largest (see Appendix C). 

2. Test Two to Test Seven 

Figures 15 to 20 show test two to test seven for increasing heat 
flux with pure R-114. Also shown for comparison in each figure is the 
Churchill/Chu [Ref. 31] correlation although this is only truly valid for 
a single tube in an infinite pool. Figure 15 shows good agreement with 
the Churchill/Chu correlation in the natural convection (NC) region and 
shows no effect of the lower tube on the upper tube. The incipient point 
occurs approximately at the same wall superheat for both tubes; once 
boiling the lower tube has the higher heat—-transfer coefficient. This is 
contrary to the results obtained by Chilman [Ref. 3] and Anderson [Ref. 
13] with pure R-113, where the higher tube had the better heat transfer. 
The reason for this difference is not known, but may be due to the 
explosive nature of incipience for R-114 compared to the more gradual 
partial incipience for R-113. For R-114 experiments, the pool was 
subcooled by 1 °C for only about hour while for the R-113 experiments, the 
pool was left in a subcooled state since the previous experiment (for R- 


113 experiments, the pool is ‘heated up’ to saturation conditions). This 
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difference in subcooling may significantly affect the nature of the 
observed incipience. Further research should be conducted at _ the 
incipient point to address some of the questions. Figure 16 for three 
tubes activated shows similar behavior as Figure 15 (ie. during nucleate 
boiling, the lowest tube has the highest heat-transfer coefficient). 
Also, Figure 16 shows that the lowest tube seems to nucleate last. 
Figures 17 and 18 for four and five tubes activated show similar trends in 
both the NC and boiling regions. It appears that the tubes nucleate in 
order down the bundle (ie. the top tube nucleates at the lowest wall 
superheat and the bottom tube nucleates at the highest wall superheat). 
For tests six and seven (Figures 19 and 20 respectively), the maximum 
controllable heat flux was less than tests one through five due to the use 
of smooth tubes in the condenser limiting the condensate rate (and hence 
pressure) in of the vapor space (if enhanced tubes had been used, this 
could have been increased). Figure 19 (test six) is consistent with the 
above trends. Furthermore, the effect of activating the whole bundle 
seems to cause the lowest tube to nucleate at a lower heat flux and wall 
superheat. It also appears that there may be some influence of tube 
position in the boiling region; however, this is probably due to 
inaccuracy in the wall temperature measurements (see Appendix C). When 
the simulation heaters are also activated (test seven, Figure 20) the 
trends in the NC region are similar (ie. no effect of lower tubes on upper 
tubes). Full nucleation of the bundle however, seems to occur earlier. 
Figure 21 shows the data from tube one for all seven tests with 
increasing heat flux. This figure is of more fundamental interest as it 


shows the same tube under different bundle conditions. It therefore gives 
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a better direct comparison of the effect of heated lower tubes as any 
uncertainty in the tube wall measurements (the largest error in the 
experimental data) is the same for each test (ie. any effects seen in the 
data are bundle effects). In the NC region, tube one alone (test one) is 
somewhat different. This may be due to ‘expansion’ of the flow as it 
leaves the top of the bundle (ie. where the velocity of the flow has 
slowed down) or due to convective effects by the addition of another tube. 
For test two, Figure 21 shows an effect of the lower tube on tube one 
performance in the NC region. For test three and all subsequent tests no 
further improvement is seen. It appears, therefore that in the natural 
convection region, an upper tube is affected by a lower tube directly 
below; however, when additional lower tubes are heated, there is no 
further increase in performance of the top tube. There is also no effect 
on the incipient point (apart from test one mentioned above). In the high 
heat flux boiling region, there is also little enhancement due to the 
lower tubes. This is to be expected in an enhanced tube bundle, where the 
total heat transfer at high heat fluxes is primarily due to nucleation 
from the tube surface itself, rather than from convection around the 
surface from the tubes below. 

Figures 22 to 27 show tests two to seven for decreasing heat flux 
in pure R-114. When comparing tube one to tube two, Figure 22 shows a 
increase in heat transfer performance of tube one by tube two in the 
boiling region while Figure 15 (increasing heat flux) showed the opposite 
effect. The most probable reason for this crossover is that these two 
experiments (TBI0002) and (TBDOO02) were conducted on different days and 


startup procedures were slightly different. For TB10002 (increasing heat 
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flux), the test was conducted as outlined in Chapter IV section F (ie. 
subcooled to 1 °C, gradually heated up with data taken over a period of 
approximately 4 hrs). For TBDOOO2 (decreasing heat flux), the pool was 
subcooled to 1 °C, and then tube one and two were turned on to the highest 
heat flux (10° W/m?) and allowed to heat up for 30 minutes. The total time 
boiling for TBIO002 at the highest heat flux was therefore less than for 
TBDOOO2. For all other experiments, increasing followed by decreasing 
runs were conducted on the same day approximately 15 minutes apart. More 
research should be conducted to investigate nucleation site 
activation/deactivation. 

At high heat fluxes for test three, Figure 23 shows no heat 
transfer performance improvement between the tubes. However, at low heat 
fluxes, there does appear to be an improvement on tube one and two from 
tube three. Figures 24 to 2/7 show similar trends at high and low heat 
fluxes (ie. the top tubes are further enhanced by lower tubes at low heat 
fluxes). It should also be noted that the lowest tube in any specific 
decreasing heat flux test had the worst performance. 

Figure 28 compares tube one only for tests one to seven for 
decreasing heat flux. As stated above, there appears to be a definite 
tube enhancement at low heat fluxes with little or no enhancement at high 
heat fluxes. This is probably due to convective effects which tend to 
increase the heat transfer performance of the upper tubes due to the 
presence of lower tubes (ie. bubbles coming from the lower tubes impinge 
and slide over the upper tubes and increase the heat transfer). At high 
heat fluxes, on the other hand, all the tubes are nucleating so vigorously 


that these ‘sliding’ bubbles have little or no noticeable affect on the 
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overall performance. This supports the hypothesis of Cornwell [Ref. 7] 
that total heat transfer in a bundle is due to a summation of convective 
and nucleation heat transfer phenomena. Similar trends were found by 
Anderson [Ref. 13] and Akcasayar [Ref. 25] for smooth and finned tube 
bundles respectively (using the same apparatus) and also by Arai et al. 
[Ref. 20] for a Thermoexcel-E tube bundle. However, Akcasayar [{Ref.] did 
not find such an enhancement effect for a High Flux tube bundle indicating 
that at low heat fluxes, a porous coated surface already has a significant 
number of active nucleation sites such that impinging bubbles from below 
have little or no added effect. Turbo-B is more similar to a Thermoexcel- 
E surface and at low heat fluxes, these two types of surface obviously 
exhibit different nucleation characteristics to those of a porous coated 


tube. 


D. R-114/0IL MIXTURES TURBO-B TUBE BUNDLE EXPERIMENTS 
1. Tests with 1% and 2% oil 

Only four experiments were conducted with a 1% and 2% R-114/oil 
mixture. These were tests one and seven for both increasing and 
decreasing heat flux; no experiments were conducted for tests two through 
six. Figure 29 shows test one at 1% oil concentration for increasing and 
decreasing heat flux, clearly showing a hysteresis pattern. Compared with 
pure R-114 (Figures 13 and 14), Figure 29 shows no apparent effect of oil 
on the heat transfer in either the NC or boiling regions. 

Figure 30 shows test seven for increasing heat flux (1% oil). 


Compared with Figure 20 (pure R-114), there are similar trends (ie. no 


49 


effect in the NC or boiling regions). The tubes again appear to be 
nucleating ‘in order’ down the bundle, as found with test seven in pure R- 
114. For decreasing heat flux, Figure 31 shows similar trends to test 
seven in pure R-114 (Figure 27). Thus a 1% oil concentration appears to 
have little or no effect on bundle performance for both increasing and 
decreasing heat flux. 

Figure 32 shows test one for a 2% oil concentration for 
increasing and decreasing heat flux. Figures 33 and 34 show test seven 
for increasing and decreasing heat flux respectfully for 2% oil 
concentration. All three graphs (Figures 32 to 34) are similar to those 
for pure and 1% oil concentrations, showing that 2% oil also has little 
effect on overall bundle heat transfer performance. 

2. Tests with 3% oil 

Nine experiments were conducted with a 3% R-114/oil mixture. In 
addition to test one and seven (conducted for both increasing and 
decreasing heat flux as before), tests two through six were conducted for 
decreasing heat flux only. Figure 35 shows test one with 3% oil for both 
increasing and decreasing heat flux. Again, the figure clearly shows 
hysteresis effects between the increasing and decreasing experiments. As 
with previous oil percentages (Figures 13, 14, 29, and 32), it shows there 
is no apparent effect of oil in the NC region, similar to previous test 
one data for other oil percentages. However, at the highest heat flux 
(100 kW/m?) there is an increase in the heat transfer of about 10%. This 
is similar to the increases found by Burkhardt and Hahne [Ref.23] in a 


finned tube bundle and Arai et al. [Ref. 20]. Figure 36 shows test seven 
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for increasing heat flux. The Churchill/Chu [Ref. 31] correlation for 
pure refrigerant is plotted for comparison only. Agreement is good, 
demonstrating that in the NC region, in addition to there being no effect 
of tube position, there is also no apparent effect of oil concentration on 
the heat-transfer coefficient. As before, the tubes appear to be 
nucleating ‘in order’ (ie. top tube nucleates first with the bottom tube 
nucleating last). 

Figures 37 to 42 show data from tests two to seven for decreasing 
heat flux only. All show no effect of lower tubes on upper tubes (within 
the bundle) in the boiling region at high heat fluxes (the small amount of 
scatter is probably due to inaccuracies in the wall temperature 
measurements). Each successive figure shows that the lowest tube has the 
lowest heat-transfer coefficient; this tube is then enhanced by the 
activation of tubes below it. Again it should be noted that the 
experimental uncertainty is larger at low heat fluxes. At all oil 
concentrations, tube five is seen to have the lowest heat transfer 
performance. According to Chilman [Ref. 3], tube five had the highest 
uncertainty in the wall temperature measurements and this might be the 
cause of this discrepancy. 

If one compares Figure 42 for 3% oil with Figure 28 for pure R- 
114, it can again be seen that there is a small increase in the bundle 
heat-transfer coefficient for the R-114/oil mixture at the highest heat 
fluxes. For all tests with oil added, significant foaming was observed at 
the pool surface and this may be the cause of this increase in heat 
transfer. Schlager et al. [ref. 21] in their review article point out 


that for certain conditions (typically low pressure and high heat flux), 
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the heat-transfer coefficient increases at low oil concentrations; they 
attributed this to foaming. Figure 43 compares test one to seven for tube 
one for decreasing heat flux. As before with pure R-114 (Figure 28), 
there appears to be a definite increase in performance of the upper tubes 
by lower tubes at the low heat fluxes due to convection effects, with 
little or no such increase at high heat fluxes. 

Ore Tests with 6% oil 


The same nine tests as with 3% oil were conducted with a 6% R- 


114/oil mixture. Figure 44 shows test one for both increasing and 
decreasing heat flux. It clearly shows a hysteresis ‘loop’ between 
increasing and decreasing experiments. In comparison with other oil 


concentrations, the point of incipience occurs at a slightly lower heat 
flux. There also appears to be a small degradation in performance (10-15% 
compared with 3% oil concentration) at the highest heat flux (100 kW/m?) 
due to the oil, but there is no apparent effect in the NC region. Figure 
45 shows test seven for increasing heat flux. As before, there is no 
apparent effect of oil on the heat transfer in the NC region and the tubes 
appear to be nucleating ‘in order’. The point of incipience also seems 
unaffected by the presence of the oil. 

Figures 46 to 51 show tests two to seven for decreasing heat 
fi, At the highest heat fluxes, there is a similar small degradation 
in the heat transfer as found with test one (10-15%) when compared with a 
3% oil concentration (Figures 37 to 42). When compared with pure R-114 
(Figures 22 to 27), there is neither enhancement nor degradation, 


indicating that any enhancement provided by 3% oil is offset by 6% oil. 
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At low heat fluxes, the data are not only very similar to that for pure 
refrigerant, but also to the other R-114/oil mixtures (ie. at low heat 
fluxes, there is no effect on heat transfer at any oil concentration). 
Figure 52 compares tests one to seven for tube one for decreasing heat 
flux. As before, there appears to be the same convective enhancement at 
low heat fluxes with no enhancement (due to the successive activation of 
lower tubes within the bundle) at high heat fluxes. 
4. Tests with 10% oil 

The same nine tests were repeated for an R-114/oil mixture with 
oe O11. Figure 53 shows test one for both increasing and decreasing 
heat flux. Incipience occurred at a slightly heat higher flux than both 3% 
and 6% oil concentrations indicating that there appears to be no 
systematic increase or decrease in this point with increase in oil 
concentration. More importantly, there is a significant decrease in the 
heat transfer at the highest heat fluxes (20%) when compared with pure R- 
114. This is probably due to the re-entrant channels becoming ‘clogged’ 
with oil as the R-114/oil mixture is ‘transported’ to the surface at a 
high rate. Figure 54 shows test seven for increasing heat flux. As 
before, this shows that the NC region is unaffected by either oil 
concentration or lower tubes in the bundle. At the highest heat fluxes 
available (40 kW/m’) there appears to be little decrease in the bundle 
performance (when compared to pure R-114) due to the oil. This indicates 
that at typical evaporator operating heat fluxes, the presence of oil does 


not significant effect the heat transfer enhancement process. At higher 
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heat fluxes, however, the effect of oil appears to be very significant as 
seen in Figure 53. 

Figures 55 to 60 show test two to test six for decreasing heat 
flux with 10% oil. At low heat fluxes, there seems to be no effect of the 
oil on the local heat transfer performance. However, at high heat fluxes, 
there is a significant decrease in performance. Interestingly, if one 
compares Figures 53, and Figures 55 to 58 at high heat fluxes, the lowest 
activated tube in the bundle is significantly degraded. The effect of 
activating a lower tube significantly enhances the heat transfer from the 
tube directly above and (to a lesser degree) the tubes even higher in the 
bundle. This may be due to the vigorous boiling action of lower tubes 
partly ‘scouring’ the oil rich layer which ‘blankets’ the upper tubes. 
This effect was also noticeable with the High Flux bundle (Akcasayar [Ref. 
25% 

Figure 6] compares tube one for tests one to seven for decreasing 
heat flux. If one compares Figures 26 (0%), 43 (3%), 52 (6%) and 61 (10%) 
for tube one for all seven tests, it is clear that at low heat fluxes, the 
heat transfer coefficient is similar, regardless of oil concentration. 
Furthermore, convective effects are consistent and provide similar 
enhancements in heat transfer performance for all concentrations. At high 
heat fluxes at 0, 3, and 6% oil concentrations, there is little 
enhancement due to activation of lower tubes. However, at 10%, there does 
appear to be a small heat transfer enhancement due to activation of lower 
tubes. This was attributed above to increased ‘scouring’ of the oil from 
the vicinity of the Turbo-B surface by the increase in bubble activity as 


more tubes are activated within the bundle. However, for a practical 
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operating heat flux range between 15 and 30 kW/m’, there is no significant 
degradation in heat transfer performance for an oil concentration of up to 


103%. 


E. COMPARISON OF R-114/0IL MIXTURE EXPERIMENTS 

Figures 62 and 63 compare tube one from test one for increasing and 
decreasing heat flux for all oil concentrations. Figure 62 shows no 
effect of oil in the NC region, but some degradation in the boiling 
region. The correlation of Churchill/Chu is included for comparison. The 
incipient point appears relatively random, indicating no early or delayed 
nucleation caused by the presence of oil. Figure 63 shows similar 
degradation with a significant effect of the oil (20% decrease in the heat 
transfer from 0% to 10% oil) at the highest heat fluxes. 

Figures 64 and 65 compare tube one from test seven for increasing 
and decreasing heat flux for all oil concentrations. As compared to tube 
one test one (Figures 62 and 63), the presence of oil has some degradation 
effect on the heat transfer performance (15%) in the NC region; this may 
be due to a change in the mixture properties which would tend to increase 
the wall superheat slightly as shown. As with Figure 62, Figure 64 shows 
that the incipient point appears relatively random. Figure 65 shows 
similar trends to Figure 63 (ie. no effect of oil at high heat fluxes), 
but shows convective effects at low heat fluxes. This was expected and 
previously reported (Figures 28, 43, 52, and 61). 

Figures 66 and 67 show the average bundle heat-transfer coefficient 
(ie. an average of all five instrumented tubes) as a function of heat flux 


for test seven for increasing and decreasing heat flux respectively at all 
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oil concentration. The data are from the same data set as that shown in 
Figures 64 and 65. Comparing Figures 64 and 66, degradation is seen in 
the NC region (15%) due to the change in mixture properties as mentioned 
above. However, at a practical operating heat flux range between 15 and 
30 kW/m?, the presence of up to 10% oil causes no degradation in bundle 
performance as seen. These trends over this heat flux range were similar 
to that found for a High Flux tube bundle (Akcasayar [Ref. 25)])). 
Comparing Figures 65 and 67, similar trends (ie. no significant effect of 
oil) are found. However, due to the limit in maximum controllable heat 
flux for test seven, data at ‘higher’ heat fluxes ( up to 10° W/m*) could 
not be obtained (as mentioned earlier). As shown in Figures 62 and 63, 
there may be a significant degradation in the heat transfer performance at 


these higher heat fluxes. 


ide COMPARISON WITH PREVIOUS NPS DATA 

Figure 68 shows a comparison between the present data for a Turbo-B 
bundle, the data of Anderson [Ref. 13] for a smooth tube bundle and the 
data of Akcasayar [{Ref. 25] for both a 19 fpi and High Flux tube bundle in 
R-114. For clarity, only test one (tube one) for a decreasing heat flux 
in pure R-114 has been shown. Figure 68 shows that the Turbo-B tube has 
a significantly lower heat-transfer coefficient than the High Flux tube at 
all heat flux. This is surprising since Sugiyama [Ref. 34] showed that in 
the single tube apparatus, the Turbo-B tube was the best performer. The 
reason for this difference in behavior is not known. Also, all of the data 
appear to be parallel to each other; one may expect the enhanced tubes to 


have a different slope to a smooth tube due to the greater amount of 
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nucleation. At the highest heat fluxes, Turbo-B and the 19 fpi tube 
appear to have a similar heat transfer performance. However, the finned 
tube heat flux is based on the root diameter; if the actual finned area 
had been used, then the heat flux would be significantly lower. Figure 68 
shows that the heat transfer enhancement given by the Turbo-B tube when 
compared to a smooth tube is about three at high heat fluxes and increases 
to about five at low heat fluxes. 

Figures 69 and 70 compare the average bundle heat-transfer coefficient 
(ie. test seven) for a given oil percentage to that with no oil for all 
four tube bundles at heat fluxes of 15 and 30 kW/m’ respectively. 15 and 
30 kW/m? were chosen as being representative of the lower and upper limits 
of heat fluxes used in practical Naval evaporators. At 15 kW/m*, Figure 
69 shows large enhancements for the smooth and finned tube bundles for all 
oil concentrations, especially at the lower oil concentrations. However, 
the Turbo-B and High Flux tube bundles show a degradation in the heat-— 
transfer coefficient at all oil concentrations (approximately 5-10% at low 
oil concentrations dropping to nearly 25% at 10% oil for the High Flux 
bundle). At the higher heat flux (30 kW/m?) , Figure /70 shows similar 
trends as Figure 69. However, the High Flux bundle now exhibits a 40% 
decrease in the average bundle heat-transfer coefficient at 10% oil. 

Figures /1 and /2 compare the average bundle heat-—transfer coefficient 
(ie. test seven) for each enhanced tube to that for the smooth tube bundle 
(tested by Anderson [Ref. 13]) for all oil concentrations at heat fluxes 
of 15 and 30 kW/m’ respectively. At 15 kW/m* with pure R-114, Figure 71 
shows an enhancement factor of 3.7 for the Turbo-B tube bundle. This 


enhancement decreases slowly with increasing oil percentage to a factor of 
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about 2.5 at 10% oil. This agree very closely with the 19 fpi bundle. 
The High Flux bundle exhibits much larger enhancements, from over 6 at 0% 
oil to just over 3 at 10% oil. At 30 kW/m? with pure R-114, Figure 72 
shows an enhancement of 3.8 for the Turbo-B tube bundle, decreasing to 
about 3 at 10% oil. This again agrees closely with the 19 fpi bundle. It 
should be noted that the High Flux bundle enhancement has decreased to a 
similar value and gets worse than the other bundles as the heat flux is 


further increased. 
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Figure 12. Performance of Test One For Preliminary Experiments 
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Figure 13. Performance of Test One at Various Tube Positions for 
Increasing Heat Flux 
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Performance of Test One at Various Tube Positions for 
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Figure 14. 
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Figure 15. Performance of Tubes 1 and 2 for Increasing Heat Flux in Pure 
R-114 
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Figure 16. Performance of Tubes 1, 2, and 3 for Increasing Heat Flux in 
Pure R-114 


63 





Sol a eee 
ee a | a 
aa ed 
Hale 
Me coo 
Mo MECC 


HELL Mt 
ANC 








TWALL - TSAT (K) 


ttt i . 
Cet e 0A 
tno 0200600? 


999990? 










oO oO oO 2 
Oo Oo ) © 
O © © = 
>) © = 

© rT 

= 


(ZIW/M) XN1d LVSH 


Figure 17. Performance of Tubes 1, 2, 3, and 4 for Increasing Heat Flux 
in Pure R-114 
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Figure 18. Performance of All Five Tubes for Increasing Heat Flux in 
Pure R-114 
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Figure 19. Performance of All Five Tubes with Active Pairs for 
Increasing Heat Flux in Pure R-114 
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Figure 20. Performance of the Bundle 
Increasing Heat Flux in Pure R-114 
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Figure 21. Comparison of Tests One to Seven for Tube 1 for Increasing 
Heat Flux in Pure R-114 
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Figure 22. Performance of Tubes 1 and 2 for Decreasing Heat Flux in Pure 
R-114 
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Figure 23. Performance of Tubes 1, 2, and 3 for Decreasing Heat Flux in 
Pure R-114 
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Figure 24. 
in Pure R-114 


4. Performance of Tubes 1, 2, 3, and 4 for Decreasing Heat Flux 
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Figure 25. Performance of All Five Tubes for Decreasing Heat Flux in 
Pure R-114 
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Figure 26. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in Pure R-114 
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Decreasing Heat Flux in Pure R- 
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Figure 28. Comparison of Tests One to Seven for Tube 1 for Decreasing 
Heat Flux in Pure R-114 
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Figure 30. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 1% Oil 
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Figure 31. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 1% Oil 
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Figure 32. Performance of Tube 1 for Increasing/D i : 
R-114 with 2% Oil ng/Decreasing Heat Flux in 
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Figure 33. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 2% Oil 
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Figure 34. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 2% Oil 
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Figure 35. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 3% Oil 
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Figure 36. Performance of the Bundle 
; with Simulati 
Increasing Heat Flux in R-114 with 3% aes on Heaters for 
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Figure 37. Performance of Tubes 1 and 2 for Decreasing Heat Flux in R- 
114 with 3% Oil 
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Figure 38. Performance of Tubes 1, 2, and 3 for Decreasing Heat Flux in 
R-114 with 3% Oil 
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Figure 39. Performance of Tubes dices cer 
in R-114 with 3% Oil 
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Figure 40. Performance of All Five Tubes for Decreasin 


114 with 3% Oil 
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Figure 41. Performance of All Five Tubes with Active Pairs 
Decreasing Heat Flux in R-114 with 3% Oil 
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Figure 42. Performance of the Bundle with Simulation Heaters 
Decreasing Heat Flux in R-114 with 3% Oil 
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Comparison of Tests One to Seven for Tube 1 for Decreasing 


Heat Flux in R-114 with 3% Oil 


Figure 43. 
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Figure 44. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 6% Oil 
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Figure 45. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 6% Oil 
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Figure 46. Performance of Tubes 1 and 2 for Decreasing Heat Flux in R- 


114 with 6% Oil 
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Figure 47. Performance of Tubes Lee 
R-114 with 6% Oil 
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Figure 48. Performance of Tubes jae”) ai 


and G 
mR -114 with 6% 011 n for Decreasing Heat Flux 
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Figure 49. Performance of All Five Tubes for Decreasin 


114 with 6% Oil 
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Figure 50. Performance of All Five Tubes with Active Pairs for 
Decreasing Heat Flux in R-114 with 6% Oil 
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Figure 51. Performance of the Bundle with Simulation Heaters for 
Decreasing Heat Flux in R-114 with 6% Oil 
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Figure 52. Comparison of Tests One to Seven for Tube 1 for Decreasing 


Heat Flux in R-114 with 6% Oil 


of 


- TSAT (K) 


TWALL 


CE CCC Ce ee a et 
CIC CC CC I, Ce 


(ZW/M) XN1d LV3SH 






100000 
0000 
1000 


Figure 53. Performance of Tube 1 for Increasing/Decreasing Heat Flux in 
R-114 with 10% Oil 
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Figure 54. Performance of the Bundle with Simulation Heaters for 
Increasing Heat Flux in R-114 with 10% Oil 
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Figure 55. Performance of Tubes 1 and 2 for Decreasin 


114 with 10% Oil 
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Figure 56. Performance of Tubes sa 


R-114 with 10% Oil 
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Figure 57. Performance of Tubes 1, 2, 3, 


in R-114 with 10% Oil 
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Figure 58. Performance of All Five Tubes for Decreasing Heat Flux in R- 
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114 with 10% Oil 
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Figure 59. 
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Figure 60. 


Decreasing Heat Flux in R-114 with 10% Oil 
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Figure 61. 
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Comparison of Test One for Increasing Heat Flux in R-114 /Oil 
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Figure 63. Comparison of Test One for Decreasing Heat Flux in R-114 /Oil 
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Figure 64. Comparison of Tests One to Seven Tube One for Increasing Heat 
Flux in R-114/0il Mixtures 
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Figure 65. Comparison of Tests One to Seven Tube One for Decreasing Heat 
Flux in R-114/0il Mixtures 
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Figure 66. Mean Bundle Heat-Transfer Coefficient for Increasing Heat 
Flux in R-114/0il Mixtures 
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Figure 67. Mean Bundle Heat 
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Figure 68. Test One Comparison of Turbo-B, Smooth, Finned, and High Flux 
Tube Bundles for Decreasing Heat Flux in Pure R-114 
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Figure 69. Ratio of Mean Bundle Heat~Transfer Coefficient With Oil See 
Heat-Transfer Coefficient Without Oil for Different Oil Percentages at 
a Heat Flux of 15 kW/m? 
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Figure 70. Ratio of Mean Bundle Heat-Transfer Coefficient With Oil to 


Heat-Transfer Coefficient Without Oil for Different Oil Percentages at 
a Heat Flux of 30 kW/m? 
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Figure 71. Ratio of Mean Bundle Heat-Transfer Coefficient of Enhanced 
Tube to Heat-Transfer Coefficient of Smooth Tube for Different Oil 
Percentages at a Heat Flux of 15 kW/n? 
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Figure 72. Ratio of Mean Bundle Heat-Transfer Coefficient of Enhanced 
Tube to Heat-—Transfer Coefficient of Smooth Tube for Different Oil 
Percentages at a Heat Flux of 30 kW/m? 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Nucleate boiling data of R-114 at atmospheric pressure were obtained 
using a small bundle of Turbo-B copper tubes. The data were obtained for 
both increasing and decreasing heat flux and at different oil 
concentrations. Based upon the results pertaining to this particular 
bundle and apparatus, the following conclusions may be made: 

I Natural Convection Region 

a. For a single upper tube, a second lower tube directly below 
when turned on does increase the heat transfer performance of the upper 
tube, however when additional lower tubes are heated no net increase in 
performance occurs. 

b. The presence of heated lower tubes in the bundle reduces the 
incipient boiling point of the upper tubes and the tubes tend to nucleate 
‘in order’ (ie. top tube first, bottom tube last). 

c. The effect of adding oil to the refrigerant (up to 10%) 
reduces the heat-transfer coefficient slightly (approximately 10-15%) due 
to changes in the fluid properties. 

2% Boiling Region 

a. For pure R-114, the presence of heated lower tubes on the top 
tube causes no enhancement at high heat fluxes (> 20 kW/m*), but at low 
heat fluxes (< 20kW/m*), there is a significant enhancement due to 


convective effects. 


120 


b. At very low heat fluxes (<2 kW/m*), the presence of oil has 
little effect on the heat transfer performance of the top tube in the 
bundle. At higher heat fluxes (> 2kW/m’), the performance is enhanced by 
10-15% at low concentrations, but is degraded up to 20% at 10% oil 
concentration at the highest heat fluxes. 

c. At typical operating heat fluxes (15-30 kW/m*), the bundle 


performance is reduced between 5-15% with oil. 


B. RECOMMENDATIONS FOR FUTURE WORK 

1. Conduct experiments with varying pool height, but keep the local 
pressure at each tube constant by simultaneously varying the vapor 
pressure above the pool. 

2. Additional experiments with R-113 and R-114 should be conducted to 
investigate explosive (R-114) and partial (R-113) incipience at the onset 
of nucleation varying the time at the incipience. 

3. Some instrumentation should be added such that the flowrates 
through the bundle can be determined. From these measurements, vapor 
quality can be determined. 

4. Metal guide plates should be manufactured and placed on each side 
between the simulation tube bundle and the tube bundle itself. This 
further channels the flow of refrigerant thru the bundle at high heat 
fluxes. 

5. Attention needs to be given to the question of refrigerant 
disposal. There are reclamation projects undertaken by most 


manufacturers; however, a method still needs to be found to remove the 


Wt 


refrigerant from the apparatus into a container suitable for such 
reclamation. 

6. A high speed camera should be used to study the nucleation process 
and circulation patterns in more detail in the bundle. Neutrally buoyant 
particles might be placed in the pool to facilitate study of circulation 


patterns within the bundle. 
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APPENDIX A: LIST OF DATA FILE 


Table 3. DATA FILE NAMES FOR TURBO-B 
TUBE BUNDLE EXPERIMENTS 
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APPENDIX B: SAMPLE CALCULATIONS 


Data set number 1 Tube 1 of experiment TBD1005 (Turbo-B tube, 
decreasing heat flux, 10% oil concentration, test 5) was used for the 
sample calculations in order to validate the program used for data 


acquisition DRP4RH. The working fluid was R-114. 


1. Test tube dimensions 
Dic = 11.60 mm 
Do = 14.15 mm 
Dj = 12.70 mm 
L = 203.2 mm 
Ly = 25.4 mm 

2. Measured Parameters 
T= 10.62 26 


T2 10594 °C 


T3 = 9.96 °C 
Da 10290 e 
T5 = 10 (04> °C 
T6 = 9.07 °C 
Didiv=s 2227 °C 
Tid2 2. 215 °G 


Aas = 3.513 V 
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Vas = 3.189 V 


3. Calculations 
The heaters power is first calculated for 
q = Vas(V)xAas(V)x60(V/V) x1 (A/V) 
Note: The multiplication factors of volts and amp sensors are 60 and 
1, respectfully. 
Therefore: 
qe (> 169 )'(3 513) (60V/V) (1IA/V) 
q = 672.19 Watts 
The tube inside wall temperature is obtained from the average of all 


six thermocouple readings. 


6 


se 


Lp 


A|H 


nel 


= 1/6(10.62 + 10.94 + 9.96 + 10.90 + 10.04 + 9.07) 
0.25 °C 
The tube outside temperature is calculated by knowing the inside wall 
temperature using Fourier’s Conduction Law. Uniform radial conduction is 


assumed. 


D 
qg{in(— 





)] 
we cc 
a 2n(k.,) (L) 


ee 
A 


where the second term on the right hand side is the Fourier conduction 


term. If we define this term as 


Ihe). 


g{1n( Oe 
5 eee 
% - SEK) 





O ee esac. 


where ky, is the thermal conductivity of copper and is calculated as 


follows 


Kk pews 40 —aullomh oem a) 


cu 


Ko, = 434 Oyen 6. 112 (263.25)) 


cu 


k, - 402.28 W/mK 





now 
672.15 [ln( iris 
Fo. 7-80 
wo wi “2m (402.28) (.2032) 
T at. OGG 1) ee 


wd 


Tyo 7 9:98°C 


The liquid saturation temperature at the top of the tube bundle is 
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lai 1 1d2 


Tsat - 
2 
meee = 222! : oma 


Tsat - 2.24°C 


In order to calculate the local saturation temperature for each tube, 
correction factors are needed to account for hydrostatic pressure 
differences between the tube locations and the liquid free surface. This 


difference is calculated by: 


AP = p(g) (ht) 


For the top tube in the bundle which is 0.124 m below the thermocouple 


measuring pool temperature. 


AP = 1523.12(9.81) (0.124) 


AP = 1852.78 Pa 


For 1852.78 Pa pressure difference, corrected saturation temperature is 
obtained by adding 0.04 °C (from standard tables for R-114) to Tsat. 


Corrected Tsat is: 


Tsat. = (2.24 + 0.04)°C 


Tsat.=- 2.28°C 


is 


Therefore, the wall superheat can be obtained by the following: 


@, - Two - Tsat, 


@, - (9.98 - 2.28) °C 


Oep= 7, 70cRC 


Now that the wall superheat is known, we need to calculate the heat 
flux and the heat-transfer coefficient. To do this, we know that the tube 
is 12 inches long and is heated in a eight inch center portion of the 
tube. The unheated lengths of the tube are a one inch and a three inch 
section on opposite ends of the tube. These unheated lengths have a fin 
effect during the heat transfer process to the evaporating refrigerant. 
In order to account for this, the following procedure was adopted for both 
one and three inch sections. Calculations are shown below for the one 
inch section. Heat transfer from the unheated end is calculated as heat 
from the base of the fin: 


Qe = [(h,) () 1K.) (Ao) 110) (canine eee, 


where 
p= ™ (Do) 
=m (.01415) m 
= .04445 m 
now 


n/4(D2, — D*;) 


& 
) 


n/4(.014152 -— .012707) 


2.0578 x 1072 me 
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The corrected length of unenhanced surface at the end was calculated as 


follows 
Poa Ly + (t/2) 
m=. 0254 4 [(0.0141)5-0.0127) /2} 
= 0.0258 m 
hp 


is the natural convection heat transfer coefficient of the fin like 
ends and was calculated by using Churchill—-Chu [Ref. 22] correlation for 


natural convection for a smooth cylinder, as modified by Pulido [Ref.27]. 


fe tGpeeeee (tanh (ne.)) | 3/6 


it ; 
h, - * [0.6+.387 SAD mee iis 
° [1+ [2222 19/16] 
4 
where 


(h,) (p) ok: 


, (k.,) (A,) 


Therefore an iterative technique was necessary to calculate hp. The 


iterative technique used was to assume h, was 190 W/m*K and continue the 


iteration until successive values are within 0.001 of each other. The 


fluid physical properties are calculated at the vapor mean film 


temperature, given by the following equation. 


VSele et Two 
L 12 — > 


r5 


2.28 + 9.98 
Prim 7 ee Pca 


LOK ip Ges Ce- 27 OF ses K 


For R-114, the physical properties are given in the program by: 


Dynamic viscosity, Tgjjq in °K 


p = exp[-4.4636 + (1011.47/T itm) ]x10° 


ll 


iv 430.927x10°° kg/m s 


Specific heat, Tritm in °K 


Cy = [0.40188 + 1.65007x109(Teitm) + 1.51494x10°5(Thitm?) — 


6. 67853x10°1°(T;; me) }x10° 


C 


> = 966.31 J/kgkK 


Density, T¢ejtm in °K 


p - 16.0184533 (36.32 + 61.146414y14? + 


17.476838y2/2 + 1.119828?) 


where 
we 21 - 12.8(rfitm | 
VSS noe 
and 
y- .333 


9 = 1512.09 = 
Mm 


Thermal conductivity of R-114, T¢itm in °C 


k = 0.071 -—(0.000261) (Testq) 


eG 


k = 6.936x10°2 W/mK 


Prandtl Number 


[ (Cp) ] /k 


6.003 


Pr 


Pr 


Thermal Expansion Coefficient 


B - - (1/p) (Ap/AT) 


Peep = 2512.395 —= 


Bey = 1521.824 “2 


B = -(1/1512.395) [(.571)/(0.2)] 


B - 1.89x1073 (1/K) 


Kinematic viscosity 


_ 430.927x107% 
1512.09 


wee 2 849X107’ m-/s 


Thermal Diffusivity 


dle) 


kK 


a (pes 


2 6.936 mies 
(1512509) Secret 


@ = 4.747 x10 m/s 


Knowing the above properties, the heat-transfer coefficient hp, 


obtained by iteration 
hh = 362.57 W/mK 


Knowing this we can calculate n 


(h,) (p) 5 


am (ena 


(362.57) (44.45x10"*) ,°° 
(402.28) (3.0578x10~°) 


ne- 36.19 


then we can obtain q¢ 


Ge> ((h,)(p) (Ky) (AD)? OP canh) ae 


Gp = [((362.57) (04445) (402-28) 4a 5057b xi Ges 


(7.70) (tanh[(36.196) (0.0258)]) 


ds as Pee Wl, W 


The corresponding results for the three inch section are 
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can be 


hy = 289.47 W/mK 


1.24 W 


qf 


Therefore, the heat transfer through the heated length of the tube is 


qs = q — qf (1 inch section) -— qf (3 inch section) 
oe (67/2.19 — 1.24 — 2.51) W 
qs = 668.35 W 


and the heat flux and the heat transfer coefficient are as follows 


= qs/As 


qs/( (7) (Do) (L)) 


q 


oe. 55)/((xz) (0.01415) (.2032)) 


7.398x10* W/m? 


and finally the heat transfer coefficient 


Ps 
Bede - Tsact) 


668.35 
See a0. (7.70) 


h - 9.608x103 W/m*K 
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APPENDIX C: UNCERTAINTY ANALYSIS 


The same data run (TBD1005) was chosen for the uncertainty analysis. 
Therefore, the measured and calculated parameters found in the sample 
calculation were used in this’ section. The unceratainty analysis 
performed was for a high heat flux, but the procedure could be performed 
at any heat flux to determine the uncertainty bands. All uncertainties 
are presented as a percentage of the calculated parameter. The 


uncertainty associated with the experimental parameters is calculated from 








the equation suggested by Kline and McClintock [Ref. 35]. For example: 
R = R(xX1,X2, «06 e @ @ @ eH) 
then 
= OR 2,7_OR 2 OR 270.5 
8R Cae re) il ee et hots set arene ] 
where 
6R = uncertainty of the desired dependant variable 
Xy, = measured variables 


6X, = uncertainty in measured variables 


The boiling heat-transfer coefficient is given by 


Ane 15act) 


where 


Ds 
qtin(— 





)] 
te 


_T.- cee 
wo W2 2 (Koy) Cr) 


ral 
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In the above equation, the second term on the right hand side is 
usually called the Fourier heat-transfer conduction term. If we define 


ears as 





ye 
qg{l1n( 5 )] 


. aio 
¢ em (k,,) (L) 


and 


6, - T,.-Tsat, 


With this notation, the uncertainty in the heat-transfer coefficient 


is obtained using the following equation. 


bh _ , 8,7, , 9Asy2, oes 8b \2,, OTsat 240.5 
im list = reel = ) ae a a Tt 


where 
q= VxI 
q= V(V)xI(V)x60(V)x1 (A/V) 


and the uncertainty is 


(ge)? (Oe “0.5 
a a 


The accuracy in the voltage and current sensors are as follows 
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5Vas - +0.05 V 6Aas - +0.025 A 


From the sample calculation section 


Vas = 52169 0 Aas = 3.513 V 


Therefore, 


ig *. SVas,2., Aas .*;0.; 
= [ ¢ vas a ) J 








2 
$@ . p¢( 0205 57%, (0-025 ) }°-5 
ey ib pel las 1S ceeds 


4 - 1.72 percent 


Calculation of the uncertainty of the surface area is as follows 


Ag = (Do) CL) 








5A, ee 
A D 


"+ (SE) 08 
5 o 
Knowing the dimensions of the tube from the manufacturer and estimated 
inaccuraccies from work shop tools and human error, the unceratinty was 
calculated. 
Dimensions 
Do = 14.15 mm L = 203.2 mm 
Inaccuracies in measurements 


§D,5= O.1 mm 6L = 0.2 mm 


Uncertainty analysis performed 
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2 2 
ee (2) «( SE) jos 
A. el je 
5A 2 2 
2 ES) eae 
A, 14.15 208 | 2 
5A 
a - 0.7135 percent 


The uncertainty calculation for the Fourier conduction term given below 


bm _ Bane, Oc 2,¢5L)2)0.5 


cu 


key was calcluated using 


ky, - 434.0-(0.112(T,,)] 


fee 434.0-(0.112(283.25) } 


cu 


Koy = 402.28 W/mK 


and its uncertainty 


Mee = ((0.112(8T,,))2)°%3 


cu 


5T,; and é6Tsat are obtained using uncertainties in the thermocouple 
readings. Average wall inside temperature T,; was obtained taking the 
average of six thermocouple readings inside the tube wall. The 


uncertainty associated with this variable is 


2 
aT o4 te (2b Tei; }°°5 


wl 
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where 6T,; for each thermocouple was obtained by taking the difference 
between the measured wall temperature and the average wall temperature. 
Using this method, it has been attempted to try and take into acount some 
the uncertainty introduced by the fabrication procedure for the tube (ie. 
air gap). For this particular heat flux the following 6T,; were found for 


tube thermocouples 1, 2, 3, 4, 5, and 6. 


0.37%, ; 0.69%, , 0129) * 0, Ong 5, om ere elena toe 
OT, = Ce) et ee ee 




















On, = 42742 ¢ 


The uncertainty level for all remaining thermocouple readings (ie. 
excluding those in the tube wall which was considered to have a higher 
uncertainty) was estimated to be +/— 0.1 °C corresponding to an emf of 
approximately 4 pV. Saturation temperature was obtained by taking the 
average of two thermocouple readings and the uncertainty in this 


temperature was calculated from the following equation. 


2 
8Tsat - Eee jo-8 
$Tsat = (2¢ 224) "]08 


6Tsat - 0.07°C 


Knowing the uncertainty in the temperatures, we can now calculate the 


uncertaines in the following: 
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woe = [((0.112(8T,,))2) 5 
ee [(0.112(283.25))2] > 


5k, - 31.724 W/mK 


Now we can calculate the uncertainty in the Fourier conduction term 





Sh © (87,7, ( PXeu ya, ( 8Ly azo. 

} [ ( a re! ita = ye) 

oh _ 2 31.724 5 O12 Aonoes 
7% gcor017 2) oom | Dera as J 
= = §.072 percent 


from the sample calculations we know that 


@e= 0.2679°C dm - 0.2679(0.08072) = 0.0216° C 


therefore: 


owe. _ bg a 6A, 2 On 2 56> 2 5Tsat 210.5 
Sl iC) Se ) Th ) “+ ( 6. a car ae jeg 


= oe) a ) ee once EERE 
3 [ (0.0172) ?2+(0.07135) mya ee ene eam 
on - 7.42 percent 


Finally the calculation of wall superheat temperature 


6, a oe = Tsat 
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) =| 0.5 








60 , pecan 07 
6, Sree 





- 3.67 percent 


Table 4 shows the results of the uncertainty analysis performed. The high 
and low heat flux correspond to the approximate values of 7.5x10‘° W/m? and 
1x10? W/m* respectively. Note that the highest uncertainty (over 50%) is 
in the wall superheat at very low heat flux (1 kW/m?). This is due to the 
very low measured value of wall superheat (0.54 °C) which can not be 
accurately measured. Thus higher uncertainty occurs at very low heat 
fluxes. However, once the wall superheat gets higher (at higher heat 
fluxes) the uncertainty in wall superheat decreases significantly (to 
about 4%) indicative of the fact that the measure wall temperature is 


relatively more accurate. 
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Table 4. UNCERTAINTY ANALYSIS RESULTS 







VARIABLE HIGH HEAT FLUX LOW HEAT FLUX 
—S 7 
jut 
9% 










QO. 
: 
Dy 
= i 
7 as 
6q/q 
§D./Do 0.707% 0.707% 
ae/L 0.098% 0.098% 
6A,/As 0.714% 0.714% 
SS D231 3 
56,/8, 0 elie EE Es 
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APPENDIX D: OPERATING PROCEDURE 


A. SYSTEM STARTUP 

1. Power to the 28 kW (8 ton) refrigeration unit is provided by the 
breakers located in the main distribution panel located in the laboratory. 
These breakers were never secured. However, if power to this panel was 
lost, then these breakers must be reset. 

2. Turn the switch on the refrigeration unit control panel, located 
in front of the refrigeration unit to the "auto" position after passing 
through "on" position. This switch is always left on, unless unit was 
taken down for long repairs. 

3. Push the start button in the control box for the recirculation 
pump. This control box is located on the bulkhead above the recirculation 
pump in the outside area adjacent to the refrigeration unit. 

4. Set the desired temperature on the roughly graduated Fahrenheit 
scale on the control panel thermostat. It requires approximately one hour 
to chill the sump to — 15 °C. The thermometer in the ethylene glycol/water 
mixture (sump) must be monitored to ensure the desired sump temperature is 
attained and maintained. Slight adjustments in the refrigeration unit 
thermostat can be expected due to the coarseness of its scale. 

5. Energize the desired pumps by switching on the breakers in the main 
distribution power panel. Once the power is energized in the main power 


panel and after ensuring the pump suction valves are open, turn on the 
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pump motors by pushing down on the arm of the appropriate breaker box for 
the pumps located on the bulkhead next to the ethylene glycol/water sump. 
The pumps are marked "auxiliary condenser" (Pump #2) and “instrumented 
tube condenser" (Pump #1), respectively in the breaker box. Flow in the 
auxiliary condensate coils can be controlled with the individual globe 
valves located at the coil penetrations on the apparatus. The auxiliary 
condenser coils will produce the fastest adjustments to the system 
pressure. 

6. Energize the heater variac(s) desired by switching on the breakers 
(Bundle, and Simulation (for test 7 only)) in the main distribution power 
panel and individual breakers for each of these in the power distribution 
box (near apparatus). 

7. After ensuring that the breakers for the heated tubes desired are 
in the "on" position, follow the experimental procedures for normal 


operation outlined in Chapter IV. 


B. SYSTEM SHUTDOWN 

1. Turn all variacs to the zero position and switch off all breakers 
in the power panels. 

2. If apparatus will not be operated for an extended period, turn the 
switch on the refrigeration control panel to the "off" position after 
passing through "on". 

3. Allow the recirculation pump to operate for at least five minutes 
after switching off the refrigeration pump unit to dissipate any back 


pressure in the system. 
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4. Turn the breakers for the pumps to the off position at the switch 


boxes, and then secure the power at the main distribution power panel. 


C. EMERGENCY SHUTDOWN 


l. Secure all power at the main distribution power panel 
2. Evacuate building 


3. Call Fire Department 
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APPENDIX E: PROGRAM DRP4RH 


FILE NAME ORP4RH 
DATE November 22, 1988 
REVISED. FEB 1992 (R. RAAS) 


BEEF 
PRINTER IS 1 
Idp=d 


PRINT USING. “4%. "Select option default is 0-.*"" 

PRINT USING s6n, 90 Tahwingedats or re-srocessing previous data” ” 
PRINT USING “6X,°"1 Plotting data on Leg-Log °”” 

PRINT USING “6X, “<c Plotting data on Linear””" 

PRINT USING °6X,""3 Purge”“* 

PROUNtUottG "6X, "4 Fixup" 

PRINT USING 6%," 5 Move™~~ 

PRINT USING “6X,""6 Comb”** 

PRINT USING “6X,°"°7 Read Plot**" 


IOP [IS A PROGRAM VARIABLE TO SELECT A SU8ROUTINE 
INPUT tIdp é 

IF Ide#@ THEN CALL Main 

IF Idp=1 THEN CALL Plot 

IF Idp=2 THEN CALL Plin 

IF Idp*3 THEN CALL Purg 

IF Idp#4 THEN CALL Fixup 

IF Ido=S THEN CALL Move 

IF Idp*6 THEN CALL Comb 

IF Ide? THEN CALL Readplot 
ENO 


SU8 Main 

ICAL=THERMOCOUPLE CALISRATION 

EOMa/Ce/ CC 7) 

OIM Emf(35),7(35),01a(6),02a(6),01a(6) ,Doal(6),Lal6),Lua(6),Keua(6),&t(19), 


Ldto(23), YVolt(2),Amp(11),Twa(S),Twt(S),Theta(S), Thetab(S),0(5),Q01(5) ,Qdgp(S) 


Las) 
114Q! 
1144! 
1148 
sya 
L156 
116@! 
1164! 
1168 
ite 
1176 
1180! 
1184 
1188 
1D 
1196 
1220 
1204 
1208 
bai2 
Le6 
Loco! 


ase. | 


2 & & 


le 


DIM Htube(S),Tn(S),Tp(6) 


THERMOCOUPLE ARRAY (C( )) INITIALIZATION 

BAtaeeOO8b0Sl 25727. 94569,-761545.9295,78025595.81 
UAi(esdee (sooo09,6.d7ogGer1l, -2.56192E+13,5.9407gE+14 
READ Cie) 


PRINT HEADER AND INITIALIZE TIME CLOCK 

PRINTER IS 7@1 

BEEP 

INPUT “ENTER MONTH, DATE AND TIME (MM:00:HH:MM:SS)”,Dates 
OUTPUT OIRECTEO TO DATA AQUISITION SYSTEM (CHP 3497€A) 
OUTPUT 709: °TO*: Dates 

OUTEUIS 709. 10° 

ENTER 709:Dates 


PRINT 

PRINT “ Month, date and time :":Dates 
PRINT 

PRINT USING “!2X,""NOTE Program name ORP4RH"“** 
SEcr 


Oi oetee veariesle FOR B1SC NUMSER SOR RECORD RESPING ONLY 
INEUTe ENTER OISK NUMSER* ,On 
PRINT USING “!16X,°"Disk number = *",22",0n 
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14ic 


BEER 
Im=Q 
INPUT “ENTER INPUT MOOE (@=34S74,1°FILE) @=OEFAULT” | Ir 


INPUT MOOE ZERO IS FROM THE OATA AQUISITION SYSTEM 
IF Im*@ THEN 


BEEP 
INPUT “GIVE A NAME FOR THE RAW OATA FILE” ,D2files 
PRINT USING “J6X,"““File mame “"“,14A" O2f ir: les 


CREATE 8SOAT FILE ON THE MASS STORAGE MEOIA 
CREATE GOAT O2filet. 62 

CREATE AN INPUT/OUTPUT LINK TO OPEN FILES 
ASSIGN 6File2 TO O2files 


CREATE OUMMY FILE UNTIL Nrun KNOWN 
DifileS="“OUMMY~ 

CREATE BOAT Difile#,62 

ASSIGN @Filel TO DOlfileé 

OUTPUT @Filel :Date$ 


CREATE A PLOT FILE 

BEEP 

INPUT "“6IVE A NAME FOR THE PLOT FILE",Pfiles 
CREATE BOAT PFfile#,32 

ASSIE6N @Plot TO Pfile$ 

BEEr 


TOTC = NUMBER (TOTAL) OF DEFECTIVE THERMOCOUPEES 
INPUT “ENTER NUMBER OF OEFECTIVE TCS (O=GEFPRQUET = tare 
{DTC = LOCATION OF DEFECTIVE THERMOCOUPLE 


IF Idtc#@ THEN 
PRINT USING “16X%,""No defective TCs exist""" 
GEESE 
PRINT USING “16X,""Oefective Thermocouples Indicated by -99.99°"° 
ENO IF 


SEEC 

DEFECTIVE THERMOCOUPLES MAY BE IN CHANNELS @0-E4 

THEMOCOUPLES ARE ENTEREO AS OEFECTIVE BY COMPUTER CHANNEL NO. 
JOTC=COUNTER IN LOOP FOR OEFECTIVE THERMGGCOUREES 


IF idte 0. Then 
FOR Jdtceee 10 Idte=) 
INPUT “ENTER OEFECTIVE TC LOCATION (SY COMPUTER CHANNEL NUMEERD 


“weatet Jdtic) 


1416 
1428 
1424 
1428 
1432 
1436 
144@ 
1444 
144€ 
1452 
1456 
1462 
1464 


BEEP 
NEXT Jdtc 
ENO IF 
PRINTER IS 70] 
OUTPUT 6Filel:Ldtc(e) 
i Ime] option (THIS OPTION ALLOWS OAT& ENTRY WITK DATA FILE? 
EL oe 


BEEP 
INPUT “GIVE THE NAME GF THE EXISTING OATA FILE” ,D2fiie$ 
PRINT USING “15x. “File name ~" 9146 7072 toes 


ASSIGN B6File2 TG O2fiies 
ENTER @6File2.Nrun 


12 


1468 
es ae 
L475 
1482 
1484 
1488 
L492 
1496 
1502 
1504 
1598 
WAS) Ae 
1S16 
1520 
1Sc4 
[S28 
iSO 
1536! 
1S4@ 
1$44 
1$48 
1SS< 
1SS6 
1562 
1564 
1Sa8 
ee 
LS76! 
1S8Q 
1584 
1588 
SiS 
1S96! 
1692 
1604! 
1608 
Pele 
1616 
1620 
1624 
1628! 
1632 
1633 
1634 


ENVER SE 1 lee eae ete: ltt Boo ,Nht, Nato Net. cocr 
SEEP 
INPUT “GIVE 4 NAME FOR PLOT FILE” .Pfiles 
CREATE SUMP less co 
ASSIGN Q9Plot TO Pfiles 
PRINT USING “16X,"“"“This data set taken on y pias oLdS 
BEEP 
FRINTER IS 1 
PRINT USING “4x,""“SELECT TU6GE TYFE""" 
PRINT USING “6X,""@ SMOOTH" "* 
PRINT USING “SX,""1 FINNEOC1LS/IN) “*"" 
PRINT USING "6%." "2 KRIGH FLUX -"* 
PRINT USING “6X,"°3 TURSO-8 "°°" 
INPUT Itt 
ENO IF 
IF Im=1 THEN GOTO 1768 
PRINTER IS 1} 


IF Im22 THEN 
PRINT USING “4X,°“Select tube type°** 


PRINT USING “6X,°" @ Smooth °°" 
PRINT USING “6X,°" 1 FINNED 19/IN (OEFAULT)**" 
PRINT USING “6X,"" 2 HIGH FLUX"** 
PRINT USING “6X,°° 3. TURBO-B*°"" 
PRINT USING “6X,"°" 4 GROWTH""* 
PRINT USING "6X,°" S 6ROWTH""” 
PRINT USING “6X,"" 6 GROWTH’"” 
ITT=TUBE TYPE 
INPUT Itt 
OUTPUT OFilel: Itt 
ENO IF 
PRINTER IS 701 
Itt=2 
PRINT USING “16X,°"*Tube Type. OO wl tt 
BEEP 
Bop =@ 


INPUT “ENTER BULK OIL % (ODEFAULT20%) *, Bop 
OUTPUT @FileliBop 
PRINT USING “16X,"°"Bulk o11%2°" 00": 800 


BEEP 
Ipozl 
INPUT "ENTER POOL HEIGHT ABOVE TOP TUBE (@*LESS THEN Scm, 12Scm OR GREATER 


(OEFAULT)",Ipo 


le s5 
1636 
1637 
ieveasyy 
1640 
1644 
1648 
bes 
1653 
16S5 
L6S6 
pos 

1657 
16$8 
1660! 


OUTPUT @Filel.Ipo 

PRINT USING “16x, “Pool herght=°",00°:Ipo 

BEEP 

NHT=NUMBER OF HEATEO TUBES 

Nit=S 

INPUT “Enter number of heated instrumented tubes(default=S5)" ,Nht 
OUTPUT OFilel:Nnrt 

PRINT USING "“16X,"°"*Number of heated instrumented tubes2"*,00".Nht 
BEEP 

Ipos=l 

Ceitinit=letden INPUT “WHICH POSITION 1S THIS SINGLE TUSE ©! 190 S DEF=L)° 


Ipos=Ipos-l 
BEEP 


Joys) 


166° Natp=Number of active Gummy pairs 

1SEE Natr=@ 

1672 INPUT “Enter number of active dummy pairs (Defauit=@)" Nate 

1676 OUTPUT @Filel:Natp 

168@ PRINT USING “16X,""“Number of active dummy pairs="",00".Nate 

1684 BEEP 

1668! 

1682! NRT=NUMBER OF ADDED HEATEO TUBES TO ENHANCE BUNDLE EFFECT 

16SE Nrt#Q : 

17@@ INPUT “Enter number of added heated tubes from simulation meaters(Defaults= 
0)° het 

1704 DUTPUT @Filel:iNrt 

17@8 PRINT USING °16X,"°"Number of added heated tubes( from simulation heaters )#=" 


“, BPO saNet 
17 See eee 
17161 


172@' CORR IS CORRECTION FOR INSTRUMENTED TU8E HEIGHT 
1724 Corr#Q 
1728 INPUT “WANT TO CORRECT TSAT FOR TUBE HEIGHT (@*#YES( DEFAULT),1=NO)°,Corr 


isc IF Corr#@ THEN PRINT USING "“16X,°"°TSAT 15 corrected instrumented heat 
ed tube height’*” 
1736 IF Corr#1 THEN PRINT USING “16X,""TSAT 1s NOT corrected for instrumen 


ted heated tube height’”” 

1740 OUTPUT @Filel:Corr 

1744 BEEP 

1748! ILQV=INPUT MODE: LIQUID, VAPOR,OR LIQUID VAPDR AVERAGE 
1752 Sligqved 

1756 INPUT “SELECT (@*LIQi default), L"VAFP,2=(LIQ+VAP)/2)" ,Ilav 
1768! 

1764! DlA=Diameter at thermocouple positions (meters) 

1768 DATA .0122,0.0098,2.0120F.@.0116,8,0,2 

1772 REAO Dla(e) 

1776 Dle#OlacItt) 

178Q! 

1784! D2=Diameter to base of fins (outside dia for smooth) (meters) 
1788 DATA .@158,@.@:25,0.0155,0.01415,0,2,2 

1782 READ Orla‘) 

1796 OreDlatltt? 

1800! 

18@4! DieInside diameter of unenhanced ends (meters ) 

1806 OATA .@152,0.0109,2@.0116,8.0127,8,2@,0 

1812 READ Dial?) 

1€16 OreDiallItt) 

182Q! 

1824! Do«Outside diameter of unenhanced ends (meters) 

1€2@ OATA .@!5675,8@.0125,0.015675,.0.0:415,0,2,2 

16232 REAO Ooa(e) 

1€36 OoxDoal Itt) 

1842! 

1844! L=lLength of enhanced surface (meters ) 

1646 DATA. .2032,.27020 1.2052 .4202..,.20 seo coe 

1652 REAO Lale) 

1e56° ExlacC iti) 

1862! 

1€64' Lu=CORRECTEO Lencth of unenhanced surface at the ends (METEFE 
18S8!'LUSLFIN + THICKNESS /2 

1Gvo 0 =6OOATA .€26!,.0254,.6264.0.0258,C,C,e 

1€76 READ Lua) 

1€8@ LuFlLuat ltt?) 

18&4! 
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Pscomeve corrected denctr c'sa inen finned litle end 
Pee eee lve CG? 

Neseee0elA {O76S..07benee oc 6.0 r6S.e.¢.c 

'9ac READ Lvace) 

19@4 LvelLva(ltt) 


19@€' Keue=Thermeal Conductivity of tube 
1912' DATA 4@1,0.€,€,€,¢@,@ 

1916! READ kKeue °@' 

192@! hKeusekcuall tt) 


1924 AePJe(Gc*r-D1*2)/4 
1928 F#PIeDo 


1952 «(FH 1 

1936 Sr#@ 

1940 Sy=@ 

1944 Sxs*Q 

1945 Sxy=#@ 

1952 Repeat: | 

1956! 

196@ IF Im=@ THEN 

1964! Otld=desired temperature of liquid 

1968! Otld=47.5 IR-113 

t353 Otld=2.2 (R-114 

1972 Ido#2 

1976 ON KEY @,15 RECOVER 1952 

1982 PRINTER IS 1 

1984 PRINT USING =2x-""SELECT OPTION ~~~ 

1988 PRINT USING “6X,°°@=TAKE DATA’"” 

SIE PRINT USING "6X,°°l"SET HEAT FLUX”"*°® 

i ste|s PRINT USING “6X,°°2=SET Tsat (DEFAULT SET FOR R-114)°°" 

2800 PRINT USING “4X,"“NOTE: KEY @ = ESCAPE"”” 

2004 | Ido=desired option 

~@88 BEEP 

2012 INPUT Ido 

2@16! 

2020 BEEP 

2024! Set default value for input 

cOrs IF Ido>l THEN iocze 

Bsa! Take cate option 

2235 IF Ido=@ THEN 2442 

2042! 

2044! LOOP TO SET HEAT FLUX (FOR TOP INSTRUMENTED TUBE) 

2048 IF Ido=l THEN 

2052 Oqdp=100002 

2056 PRINT USING °4X,°°Qdp QDPsim Nrt Qdpaux 
Ceote a: 

2062 PRINT USING "4X,""C(W/m"2) (Wim?) (W/m*s) 
CW 

2064 Erre] 

2066! Reset,read channel 25-3@,automatic scaling 

fal GCS Cnanmne!] 2S=aur» amos,2Z6"51m amps, ,cl#inst volts,2@=#sim volts,29=au-, 
volts,30-34=:1nst amps 

2076 OUTPUT 7809: "4&R AFZS AL354 VRS" 

2082 FORT £10 1071) 

c@es OUTPUT 705."AS SA” 

2OBE ENTER 709. Amp! 1) 

~@92 NEYT I 

2036 FOF T=? TC 2 

-102 OUTPU. 729: °AS $4" 

2104 ENVEF .205.. Vol) 

~1@8€ NEXT 2 
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elCs FOR i=€ Tera 

Ske OUTEUT 7ee. "AS Se” 

BOS UE ENTEF 7@5.Amp( 1) 

Pasa as NEXT I 

2120 Calculate actual heat flu> 

efes OC @ =6G9Volt(@ eAmp(Ipos) 

Zice Qan(€)=Q(0)/(PIeO2eL ) 

SLoe OsimeB@elQ@eVolt() deAmp (11) 

£135 OdpsimeGsim/(FleO2*.2@32¢5 ) 

2142 Oeux*B@elBeVolti zc deAmp( 10) 

2144 OQdpaur *Qaux/(Ple.@160¢.1776e4) 

eae Qrot=Q(8)eNnt+Osim+Gaux 

21S2 Nrt=Odpsim/Qdp(Q) 

2156 IF ABS(Agdp-Ogo0p)>Err THEN 

2162 IF Agdp>Dadp THEN 

2164 BEEF 400@2,.2 

2168 ELSE 

2lfzZ SEEP 25@,.2 

2176 END IF 

218@ IF Nrt<.1 THEN Nrt=@ 

2184 IF Qdpaux<1@@ THEN Qdpaux=@ 

2188 IF Qdpsim<100 THEN Qdpsim=@ 

fice PRINT USING “4X ,2¢M2.30E, 2X) 2x, (M00. DO). 2x 2 eee 2k oes 

(@),Qdpsim,Nrt ,Qdpau> ,Otot 

216 WAIT 2 

2202 6O0TO 28076 

2204 ENO IF 

2208 END IF 

Bead 

c2lb! LOOP 70 SET Tsat 

eae IF Ido*2 THEN 

e2a4 IF Ikdt#1 THEN 2240 

2228 Geer 

Beoo! INPUT “ENTER DESIRED Tsat (DEFAULT#<47.5 C = R=1] 30m, 021d 

22o0 INPUT “ENTER DESIRED Tsat (DEFAULT=¢.2 C - R-114)" ,Otld 

2236 Ikdt=1 

2240 Cidl=@ 

ead Oid2=@ 

2246 Nn=] 

2282 Nrs#*Nr MOD 15 

e256 Nn=Nn+1 

2260 IF Nes#] THEN 

2264 PRINT USING °4X,°° ODTsat Tle! Waid2Z Tlbb 
Tvat Tvab Tlav epi 

Z2be ENO IF 

Ee FU Reac thermocouple voltages for vapor, liquid 

2. 7b OUTPUT 7@5.°AF AF@ ALS VRS" 

cao! Sample each tnermocouple 2@ times and report temp for each the 

rmocouple, vapor*@,1,¢2; IJlaiquid#584 

wend FOR :#@ TO ¢ 

mabe Suma 

2292 OUTPUT 705: "AS SA” 

2296 FOR Ji=1 TO 20 

27202 ENTEF 709.Elig 

2304 Sum=Sumtli19q 

~ 3@6 NEXT Ja 

Ba ve Emfiy :=Sum/2@ 

2216 TCT ds FNTvsviEmf(] )3 

aoe NEYT I 

meee! Compute average temperature cof liquid 
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eoe Lieve t 1s fe T Ge) )2../5 


Eee ag Comoute average temperature of vapor 
Boo Teas t=t1 CG 4101 72 

Boo TvavesT(l) 

ae Teevet (UC) ttt (eres, = 

nae [FACS bev-Ovle) 22.) THEN 

2344 TE Tas Pile THEN 

7348 BEEP 4002,.2 

239. PE SE 

egigls BEEP 2S2,.c 

< 562 ENO IF 

2364 ELSE 

Bobe IF ABS(Tlav-Otld)>.1 THEN 
Ler IF At1ld>Otld THEN 

2276 BEEF esOO0..2 

2382 ELSE 

2384 BEEP 800,.2 

ZouG ENG ie 

2392 ENO IF 

2396 END JF 

2400 Errl=Tlev-Dldl 

2404 OldleTlav 

2408 Err2=Tvav-Dld2 

Pelz Did2"Tvav 

2416 PeaNioWelNoe «44.7 CMDU0700,5%) sOtld,1(2).9¢46),1(5),Tvavi,Tvave,Tla 
Vv 

2420 WAIT 2 

2424 GOTO 2252 

242 END IF 

2432! 

2436! TAKE DATA IF Im=#@ LOOP 

2442 IF ITkole=l THEN 2452 

2444 SEEF 

2448 Ikol#=1 

24$2 DUTPUT 729; "AR AF@ ALS VRS" 
2456 FOR J=?310 5 

2462 OUTPUT 70S:"AS SA” 

2464 Sum=@ 

2468 FOR J1=i7 9029 

7a a ENTER 7@S;E 

2476 Sum=Sum+e£ 

2482 IF ]>2 THEN Et(Ji-lieE 
2484 NEXT Ji 

2488 Kdl=@ 

2492 IF I>2 THEN 

2496 EaveeSum/20 

2-502 Sum=@2. 

2504 FOF Jre@ TO 19 

2508 TF ABSC(C EMC J) J-Eave )<S.0E-E THEN 
Zoe SumeSum+Ee( Jb) 
Veils aisle 

2522 Kdlekd}+] 

2524 ENO IF 

Goac NEXT Jb 

Sooc LE Pee eee her RING USING “Oh. Reale °° OR” ek) 
Be oo US iis Oe Gass 

2542 Ecer 

S544 Baer 

~54e PFINT USING °4¥,""Too mucn scattering in data - re 


peat cete se. ~ 


Mei! 


~SS2. 
2556 
256e@ 
2-564 
<S6E 
Resign 
2576 
-58e 
Boos 
A Sse 
eae 
2596 
2622 
2604 
Z6QE 
Bead cd 
2616! 
2622 
2624 
2628 
2652 
2636 
2640! 
26441 
2648 
2652 
2656 
2662 
2664 
2668! 
2672 
2676 
2680 
2684 
Z68e 
2632 
eocu 
2702 
2704! 
2728 | 
2tie 
2716 
2722 
2724 
2728 
EToe 
27356 
2742 
2744) 
2748 
2749 
2752 
ae) 
27Ez 
154 
27E5 
2756 
ay LV 
2755S 
~7E2 


GOTO 1S8¢ 


ENC IF 
END JF 
Emf iy s=Sum/,y2Q-be] } 


NEXT ] 
OUTPUT 7@9."AR AF4C ALES YRS“ 
FOR 1=6) 1035 

OUTPUT 7@9."AS SA” 

Sums? 

FOR Jie eee 
ENTEF 7@095,.E 
Sum=Sum+& 

NEXT Ja 

Emfi lt )=Sum/S 

NEXT ] 


REAC VOLTAGES (27=Inst,2&=Sim,29=Au> ) 
OUTPUT 7@9.°AR AFL? ALZS VRS” 
FOR I=¢€ 10 2 
OUTPUT 709:"AS SA” 
ENTER 709; Ueber] ) 
NEXT I 


READ CURRENTS (3@-34=Inst tubes: 35-39sACTIVE Dummy ) 
OUTPUT 729; "AR AFZ2@ AL3S YRS“ 
FOR I=@ TO 9 
OUTPUT 70@9.°"AS SA” 
ENTER 709. Amp(1 > 
NEXT I 
Read Currents(25=Aux amps,c6=Sim amps) 
OUTPUT 709:;"AR AFZS ALZ6 VRS" 
FOR [#10 TO 11 
OUTPUT 7@9:"AS SA” 
ENTER 7@9:Amp( 1) 
NEXT I 
ELSE 
ENTER @FileciEmf(*) Volt) Amol.) 
ENC IF 


CONVERT EMF’S TO TEMP,VOLT , CURRENT 
FOR I#*@ TO 35 
TCT @FNTvsv(Emf(l)) 
IF 155 AND Idtc>@ THEN 
FOR I1*@ TO Idtc-l 
IF Lotet [i = 1-439 THEN eee 35,49 


NEXT IT, 
ENO IF 
NEXT ] 
Ntc=nr of thermocouples 
Ntc=6 


IF Ipos’@ THEN 

CCI pos © 52*V0lt(@)*Amp( Ipos ) 

TwallIpos =2 

J322 

Ndtc#@ 

FOF T=] TO Ntc 
Nn=Tpos*6+5+J) 
Jj=J3+) 
Te ABS(T(Nn): 9S THEN 

TON =-GE, 95 
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~761 
ar 6c 
S755 
=e: 
2765 
2766 
B76? 
2768 
2765 
2770 
2771! 
Pare 
2773 
2774 
2776! 


ermocouple in 


~780 
2784 
2788 
erd2 
“fests 
2820 
2804 
e828 
Z2el2 
2816 
2820 
eee 
2824! 
2828! 
2e.9! 
263< 
2836 
2840! 
28441 
2848 
eoa2 
c@56 
2860 
2864 
2868 
e872 
2876! 


3 


EeT7} 
ee 
2882 
2884! 
2886! 
we 
~896 
2922 
29Q4 
-9@6 
eae 
2516 
2520 
oes 
2928! 


eos. | 


Ndote=Ndves! 
ELSE 
Twe(Ipos )=Twa'i Ipos )+T( Nn) 
ENG) IF 
NEXT ] 
Twa(lpos )=Twa(Ipos)/(6-Nodte) 
GOTC 782e 
IF 
I,1=@ TC 4 
Q: li *BQeVol~(Gd*Amp( 11) 


END 
FOF. 


TwatAverage temperature of the wall 
Twa( li )=Q 
Ndtc#=@ 
FOF I=1 TO Nte 
No ts ecoumter lmetenpvactray. start at © (this is the first th 


the tube tant ) 
Nn#=J]1¢6+] +5 
IF ABS(T(Nn))>99 THEN 
T(Nn )#-99.99 
Ndtc=Ndtc+! 
ELSE 
Twa( li Je Twa( Ji )+TCNn) 
ENO IF 
NEXT I 
Twa( di )#=Twa(J1)/(6-Nodte ) 
NEXT Ji 
Tlav=(T(3)4+T(4))/2 
TvaveT(2) 
Tvave(7T(Q@)4T(1 047T(2))/5 


Tlav2T(S) 
Tou=Twa(Q) 
Kou=FNKeu( Teu) ITHERMAL CONOUCTIVITY OF COPRFER 
'TF CURVE FIT NOT AVAIL USE ARRAY KCUCe ) 
FOUPIER CONOUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
FOR I=@ TO 4 
Twt tl =Twatl )-Q¢] )*LOG(02/01 )/( 2ePT eK curl ) 
IF Ilqv=@ THEN Texs#Tlav 
IF Ilqv=1 THEN Texs#Tvav 
IF JTIlqvecl THEN Texs#(Tlav+T(2))#.5 
IF Corr#=] THEN Thetabd(] )=Tw(] )-Texs 
IF Corr=@ THEN Thetabd(] )=Tw( J] )-(Texs+.@S56+1¢.129) 'tR-114 
IF Corr=@ AND Ipo=l THEN Tnetab(I )=TwlI -(Texs+.Q@54+1]¢.144) '!R-]1 


IF Corr#@ AND Ipo#@ THEN Thetab(I )=Twl I )-€ Texs-1.@78+.147¢1) l= 


NEXT 1 


COMPUTE SVARITOUS=PFGPERTIES 
Ttvim=t)aCletexs)e.5 IFILM TEMPERATUPE 


Rho=FNRho(Tfilm) 
MuEFNMu( Tfilm) 
beCNt(Tfi lm) 
CpzFNCp(Tfilm) 
Beta=FNBeta(Tfilm) 
Ni=Mu/Fhe 
Alphssk/(PheoeCp } 
Fr=Ni/Alphe 


IDENSITY 

IN TSEOSI TY 

PIREP HAS CONDUCT LU TIY 
PeeGe (PPC Real 
PTRERMAL. EYPANSION 
DUNE Me Tie Vi SCC Si. Ty 
DIG E MACS EIEEUEL UT TY 
1PRANOTL 


COMPUTE NATUPAL-CONVETTIVE HEAT-TPANSFER COEFFICIENT 


Isis, 


cose, FOF UNENHANCEC END(S) 


-942 LuScwauntt! 

2944 Hbar=190 

-948 Fee(HpareFP/(kouen))’ .SeLy 

Soo Tann=FNianh( Fe) 

2956 Theta(lIpos )=Tnetab(Ipos )eTann’/Fe 

“960 Xxe(9,61¢BetaeThetab!Ipos ec "seTanh/(FeeNieAlpma)> . 166667 

2864 Yye(1l+( .§59/Pr 5709416 )9°( 8/27) 

-36E HbarceK/Doel .6+.,3687¢exx/YVy >*s 

ye IF ABS( (Hber-Hberc )/Hbarc )>.@@) THEN 

2376 Hbar®=(HbartHbarc )e.S 

- 982 GOTC 2aee 

~984 ENO IF 

2986! 

Fay || COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S ) 

299g Q1(@ d#(Thetab(Ipos )«Tanh del (HbareFekKcueh)”.5) 

2282 Gq=01(8)+0q 

3004 Z=i+1 

3Q@2E IF 2Z=*1 THEN 

3@12 Lu#lLv 

3@16 GOTO 2944 

3822 END IF 

302 Z=2 

3028 Qlpct=Qq/O(lIpos)? 

3222 Qq=2 

3@36 As=P] «D2 L 

3242 FOR 11#@ TO 4 

3044 OlCliv=6lpetsOCli) 

3048 Qdpo( 11 #(Q¢ 11 )-O10131))/As 

3052 Htube(I1 )©Qdp(I1)/Thetab(11 ) 

3Q56 NEXT Ti 

3262 PRINTER IS 781 

3064! 

3Q@68! RECORD TIME OF DATA TARING 

3072 IF Im=@ THEN 

3@76 OUTPUT 7809; "TO" 

3282 ENTER 789: Tolds 

3204 ENG IF 

3@eS |! CHURCHILL/CHU CORRELATION FOR NATURAL CONVECTION REGION 

2087 Raea=S.8leBetarThetab(Ipos )#(00)°3¢Rno/( MueAlpna ) 

3e@88 Denom=(1.4(.559/Fr )*(9/16))*( 16/9) 

3@89 Nuch#=(.6+.3874#(Raa/Denom)*(1/6))°2 

3092 Qch=*KeNuch*ThetabllIpos )/(02) 

3Q@91 | 

3@92! OUTPUT DATA TO PRINTER 

3@96 PRINTEF I5 72} 

2100 PRINT 

2124 PRINT USINS “10X,""Date Set Number = °",000,2X%,14A":J,Told$ 

218 PFINT 

Bile PRINT USING “10x,"" Tvl lve Tee Ital VleZ elas 

Tvav Tice: s °° 

eile PRINT USING "10% ,8(MD0.00,2% 2°: 10@)-101), 1(2), Mstislav as, ler 

S520 PRINT 

2124 PRINT USING “6%,°"Tube Wall Temperatures (Ges ©) Trave Odo 
H Tnhetar"*” 

2128 PPINT USING "6X%,°"% 1 Ee 5) 4 5 6B ‘Deg C) (W/r" 

~.) Vie es SES} Ges 

oes IF Ipoe>@ THEN 

3132 J = 

212e FOR Ji:*@ TO §S 
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tes Ja =Iipostes, 146) 

eae) 

NEXT Ji 

Tnn=lpos+1 

Sec LNG berm la (MOG. OO lhe a\ OOF i in al (MOU mea 


ob nite, 1 Go 


mee ip a), tpes ), 1p 4), ToS), lwat loos), Odp(lpos -.Htube( lpos +, Thetad:ipos) 


maemo.) ),lpc2),1p(5),1p(4),1p(5),Twalli),Gde(li),Htube( 11), Thetab(l2) 


See ey li MZ. 


jie: I=] 


SOIC 2177 
ENG “JF 
J3=2 
FOR I1*@ TO Nnt-] 
PORsJi=¢ 10 5 
Tet Ja d=1011 "S70 346? 


Adlai ceo 
NEXT Ja 
a = i 


FOR Jiz@ TQ 4 
Tn( Ji deleJdi 
NEXT Ja 


PRINT CUSINGm ON Oqiha a MOO. OD J eik vet Me. ODE } yin eh ( MOO OO sind 1s eT 


NEXT I 
PRINT 


PRINT USING °6X,°° Heat Flury and Tdel from Churchill/Chu Correlation 1 


ZOE),2X,1(M00.00)°: Och, Thetab(Ipos) 
PRINT 
Or =] 
IF Im*@ THEN 
BEEP 
INPUT “OK TO STORE THIS DATA SET (1#Y¥i default ),@=N)?° Oh 
END IF 
J=ethe counter for data sets 
IF Ok=1 OR Im=] THEN JeJ+] 
IF OF=1 AND Im=@ THEN OUTPUT @FileliEmf(#),Volt(*),Ampo(s) 
IF Im=] OR Di=1] THEN OUTPUT @Plot.Gdpl*),Htube(*),Thetab(e ) 
Go_on®] 
IF Im=@ THEN 
Ceee 


INPUT “WILL THERE BE ANOTHER RUN (1#Y¥( default», Q=N)?" ,6c_ 


NruneJ 
IF Gelon=€ THEN 3272 
IF Go_on<>@ THEN Repeat 
sllsle 
IF J<Nrun+! THEN Repeat 
ENO IF 
St=! 
Gee 


on 


INPUT GARE YOU SURE YOUF READY TO TERMINATE (1=Y(OEFAULT),@=NO)?" ,St 


Go_on®] 
IF Sto@ THEN 3288 
IF St|\=@ THEN GOTO 3248 
IF Im=@ THEN 
BEEP 
PRINT 


Broun teuolNs  1Oy. 7" NOVecee (22, "°° ‘deearruns were tered ar f: le 


eet las 


ASSIGN @Filel TO e 

OUTPUT @Filel: Neun-1 

ASSIGN @6File! TO Dl files 

ERNmeeweo Tel ivenet Lotetce) Itt feo. Nat ,Nete.Net. coca 
CUTRUT e's les, Cstet,Ldtet*),ltt.Boe Nht.Neto, Net, cor~ 
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Soke FOF [el iC Nurse 


Z32 ENTEF Cirle! Emite7 ve. t: * 7 ames 

Joes OUT EUs CE lel Smite geht 6. Samat «8 

goa NEXT I 

RAH re ASSIGN @File!l TO » 

Jo0 PURGE “OUMMY" 

3342 ENC IF 

2344 Cec 

3340 PRINT 

2352 PRINT USING “10,.° "NOTE ~ ""V22.°° X-¥ pairs were steceduin cic. catau: 
lle "“,1OA": J-1,Pfslef 

3356 ASSIGN @File2 TO » 

2360 ASSIGN @Plot TO » 

3364 BEEP 

3368 SUBENO 

35721 
ZA37EIXXXXXXXXXXXXXXXXXXXKXXXXXXXKXXXXXKXKXKXKXKXKXKXXKXXXKXXXXXXXXXXXXKXKXKXXKKXKXKXXXKXKKXXKXKXXXK 
2380! 


3384 DEF FNKculTcu?) 
2388! OFHC COPPER 


3392 (heleute7o > he WEL I 

3396 Kceue434-.112°Tk 125@=300K USE FOR R=] 4002.2 6 
3400! Kceu#4352.0-.1*Tk 1200-4008 USE FOR R-1e..O97°5 ¢ 
3404 RETURN Keu 

3408 FNENO 

3412! 


3416 OEF FNMu(T) 
34201 CURVE Fil OF UiscOsriy 


3424 TkeT+273.15 iC  hOok 

3426 MusEXP(-4.4636+°1011.47/Tt.3)001.0E-3 IR-114 170-360 k 
3432! Mu=.00001349(10°(SO3/(TkK-2.15))) IR113 

3436 RETURN Mu 

3442 FNENO 

3444) 


3448 DEF FNCp(T) 

o4S2) CURVE Fl OFeee 

2456 Tk=14272.15 ERNST Ls 

3460 Co=.40188+1 6S5007E-S*Tl +1 SI4S4E-6*1) "2-6 6 78SS3E—-105 1) so ihe 4 ise -40c | k 
3464! Co=(925+1.0@5°T de. O2i fR=1)3 

3468 Cp=Cp+1000 

3472 RETURN Cp 


3476 FNENO 

2480! 

3484 DEF FNRho(T) 

S48Gn weet co 7315 iC eC ake 


3492 x=]-(].8°Tk/7535.95) J TC FR 

3496 hee356.32+61.1464149x°(1/3 )416.4180159X+17.476858exX* S41 1198289x*2 
32500 RoeRo/.062426 iR=-114 

2504! Ro=1.6207479E+35-Te(2.2186546+T°2.357829)£-3) 1R-112 

250& RETURN Ro 

3512 FNENO 


2516! 

2520 DEF FNPr(T) 16000 FOF F-11)4/R-112 
S24) PreFNCpl(T)eFNMu! T)/ENE CTD 

2526) RETURN Fe 

255° FNENG 

ZSi6! 

354@ OFF FNK(T) 

3544) T<Z6@ ¢ WITH T IN C 

254E Ke. O71-.00CLE1leT 
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cisisye 
~SS6 
3560 
2564 
=56e 
BST. 
3576 
27582 
3584 
3588 
2592 
2596 
326202 
3604 
3608 
Bol2 
3616 
3620 
3624! 
3626 
S652 
3636 
3642 
3644 
3648 
3652 
$656 
3662 
3664 
3668 
B57. 
3676 
3680 
3684 
3688 
2692 
3696! 


RETURN + 
FNEND 


DEF FNianht re? 
Pel P(Fe) 
Q=EXP(-Fe ) 
Tenet > -=0)/(F+Q> 
RETURN Tanh 
FNEND 


OEF FNTvsviV) 
COM <7Ge7-CC7) 
T=C(@) 

FOR Ie] TO 7 
T=T4+C( I] )eV*] 
NEXT I 

RETURN T 
FNENO 


DEF FNBeta(T) 
Rop=FNRho( T+. 1) 
Rom=FNRmo( T-.1 ) 
Beta=-2/(Rop+Rom )#(Rop-Rom)/.2 
RETURN Beta 
FNENO 
DEF FNPoly(X) 
Comey Colv, AC10.10),6110),804),Nop,lorcnt,0p00,llog 
X] =X 
Foly=8(@) 
FOR I=] TO Nop 
IF Ilog=1 THEN X1l#LOG(X) 
Poly*Poly+B( I )#x1°I 
NEXT I] 
IF Llog=l1 THEN Poly#EXP(Poly) 
RETURN Poly 
FNEND 


S7QQIXXXXXXXXXXXKXXKKXKKKKKXKKKKXXXKX KKK XK KK KKKK KKK KK KKK KKK K KKK KKK KKK KKK KAXKKKXX 


27@4 ! 
378e 
She 
2alé 
3722 
oe 
Bice 


Res Be Aon) 
aw five 


2736 
2742 
3744 
2748 
2752 
3756 
376@ 
3764 
3768 


fa ag a 


Bio 
Z7Be 
3784 
ZTE 


SUE Poly(OfileS(e¢),No,Itn) 
Bly R010 ),5(10),Sy1(12),S4( 12), X41 100), Yy(100),xy(17) 
ote Goly7, Av1e@,1@),Ci10),804).N, lornt,050,lloo 
COM /Xxyy/ Xp(S),Yo(5) 
FOR I=@ TO 4 
B( 1 )=«@ 
NEXT I 
Ime] 
BEEF 
INPUT “ENTER OATA FILE NAME” ,Dfiles(@) 
BEEP 
INPUT "ENTER NUMBER OF ¥-Y PAIRS“ Np 
SEEP 
INPUT “LIKE TO EXCLUOE OATA PAIRS (J) #1, @=N( DEFAULT) )?” led 
IF Ted=1 THEN 
BEEP 
INPUT "ENTEF NUMSER OF FAIRE TC EE EXCLUDEO”, Ipe> 
ENO IF 
Gaston 6rile 10 Ofiles(Q) 
N=? 
BEER 
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5792 
7796 
3602 
3804 
ZEQE 
5612 
SELE 
3622 
3624 
Sece 
5652 
3E36 
3842 
3844 
3848 
3e52 
3656 
3862 
3864 
3868 
3872 
3876 
3882 
3884 
3888 
3892 
3896 
3922 
2904 
2908 
3912 
3916 
3922 
3924 
3928 
2936 
3942 
3944 
3946 
3952 
2956 
3962 
3964 
3966 


s6>> 
ww! & 


3976 
3982 
3984 
398E 
gd92 
3995 
4220 
£004 
400E 
4012 
4Q1€ 
4020 
4004 
4G2e 


INPUT "“ENTEF THE OROEF OF PCLYNOMIAL 
FOF I=@ TON 
Sy( I )=#@ 
Sx( J] #2 
NEXT I] 
IF ITed=1 ANO Im=#!} THEN 
FOR I#! TC Ipe> 
ENTEF @File.Xyle) 
NEXT I 
ENO IF 
FOF I=] TO No-Ipex 
ENTER @File:Xy(e) 
IF Opo#@ THEN 
Yexy(Itn-1)? 
KeXy(ll+Itn) 
ENO IF 
IF Opo#=l THEN 
y=Xy(StItn) 
KeXy(ll¢Itn) 
END IF 
IF Opo#2 THEN 
Yexy(S+Itn) 
Kekvcrtn=|] ) 
END IF 
IF Ilog=#=! THEN 
X*#LO6(X ) 
Y*eLOGC(Y) 
END IF 
Xx (] )=Xx 
Yy(1)#¥ 
R(Q)#Y 
Sy(@d=Sy(O)+Y 
S( 1 )=X 
Sy C1 89x01] +X 
FOR J=1 TON 
R( J d#R6J-1 )*X 
Sy( J #Sy( J +R] ) 
NEXT J 
FOR J#2 TO Neo 
Sel 7=5( 1-1 ex 
Sx( J )eSx(J 4505) 
NEXT J 
NEXT I 
Sy (Q)=Np 
FOR I#@ TO N 
CiljeSy( 1) 
FOR J=@ TON 
ACI, J d=#S-(I+J) 
NEXT J 
NEXT ] 
FOR I=@ TO N-] 
CALL Divide(]) 
CALL Subtract(!I+!} 
NEXT I 
BOND=ECINI/ACN,N) 
FOR J=@ TO N-!} 
BiN-l-JT v8C(N-J-1) 
FOF J=@ TO ] 
BUN-)-l3=-6°N-!- | - ah -1 hd 
NEXT J 
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COEFAULT=2 "NM 


)eBON- J) 


4052 BCN-J-ij mG (N-J-] > AVN-1-3 |N-J]-] 
6056 NEYT I 
4040 IF Jornt«=C THEN 
6044 PRIN? USING "12>." "EXPONENT COEFFICIENT*”" 
404 FOF J=@ TON 
4@S72 PRIN USING “J5.-00 5) .M0.70E°.1,801) 
40S6 NEXT I 
4262 PRINT ” ° 
4064 PrN ievoinG lex, OAtA POINT > } YCreccueel eo > Oi 
SCREPANCY**” 
4066 FOR I=} TO Np 
4&Q@7e Yco=B(@) 
4@7E FOF Je! TON 
42B2 Ye=Y¥c+B( J exe (19°) 
4QB4 NEXT J 
4@€6 DeYyil)-Ye 
4232 PRENT WSING ~15xX,30-4x% ~6(Mi0.SOE.1%)° . 1, Xe01),¥y( 1? y¢6.0 
4296 NEXT I] 
41@@ END IF 
41@4 ASSIGN OFile TO * 
41@8 SUBEND 
eli2! 
4116 SUE Divide(mM) 
fee 6 COM 6/Cply/ AC1@,10),C(10),8(4),N,lIornt ,Opo,llog 
4124 FOR IJ=mM TO N 
412B Ao=A(1,M) 
4132 FOR J=M TO N 
4136 ACI, J #AC1,J9/Ao 
4142 NEXT J 
4144 GC] y=C(1)7 fis 
4148 NEXT ] 
Sis2 SUBEND 
4156! 
416@ SVE Subtract(hk) 
4164 COM /Coly/ AC1@,1@),C0(10),6(4),N,Ipent ,Qpo,Ilog 
4168 FOR I=K TON 
Oe FOR J=k-1 TON 
4176 ACT JosA(K-1,J)-ACI 2) 
4\1B2 NEXT J 
4184 CCT eC h-1)-CC(]) 
4186 NEXT I 
4132 SUBEND 
4196! 
420@ SUE PlIin 
S204) GCOM /Coly/ A(10,1@),C(10),6(46),N,.lornt ,Opo,Iloee 
wees COM /Xvyy/ Xx(5),Yy(S) 
46212 PRINTER IS 725 
4216 BEEP 
ao INPUT “SELECT (@=h/hO% same tube, L=nCHFE )/h( sm)" Irt 
4224 BEEF 
GeecmeeiNe UT WHICH Tsat (1=6.7,0=-<.2)" Isat 
4232 Xmine2 
4236 Xmax=1Q 
£240 Xstep=2 
6044 IF Irtz=@ THEN 
2246 Ymin=@ 
re Ymarvel.é 
20S s Ystept.cz 
£76@ ELSE 
£io¢ Ymine@ 
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426¢ Ymar21S 


ara Veteo=s 
2276) ENG IF 
£280)" BEEF 


4284 PRINT "IN; 591.IF 2306,2200,8350C ,680C. * 
4°62 PRINT “SC @,10¢.C. 10C 1) 270. 
BOSD Sf>81Q00/(Xmar-Xmin) 

4296 Sfy=#1@@/'Ymar-Ymin? 

420@ PPINT "PU @,@ PD” 

£3@4 FOR xaeXmin TO Xmax STEP step 
4308 Ke(Xe-Xmin)*Sfx 

43i2 PRINTER ke On KT 
43216 NEXT Xe 

432@ PRINT "PA 10@,0:FU;" 

4324 PRINT “PU PA @,@ PD” 

4326 FOR YatYmin TO Ymae STEP Ystep 
4332 yo(Ya-Ymin)#Sfy 

4336 PRINT "PA @,°~:Y,°YT* 

4340 NEXT Ya 

4344 PRINT “PA @,100 TL @ 2” 

4348 FOR Xe=Xmin TO Xmar STEP Xstep 
4352 Xe(Xa-Xmin)#Sfx 

4356 PRINT "“PA°:X,",100. XT" 
436@ NEXT Xa 

4364 PRINT “PA 108,100 PU PA 100,@ PD" 
4368 FOR Yae#Ymin TO Ymar STEP Ystep 
4372 Ye(Ya-Ymin)eSfy 

4376 PRINT “PD PA 100,° .Y,°YT” 
4380 NEXT Ya 

4384 PRINT “PA 102,100 Pu" 

4388 PRINT “PA @,-2 SR 1.5,2° 

4392 FOR Xa=Xmin TO Xmax STEP Xstep 


4336 X=(Xa-Xmin )eSfy 
442Q PRINT 7FA> 3). 7.03 
4404 PRINT “CP -2,71,LB8°: Xe; °" 


4408 NEXT Xe 
44:2 PRINT “PU PA @,2° 
4415 FOR YasYmin TO Ymax STEP Yatep 


4420 IF AESCYad<].E-S THEN YesO 
b424 Y=(Ye-YmindeSfy 

4428 PRIN) Pe @ocmi 

4432 PRINT| Ce e-6 (-. 250 Ge aa 


4436 NEXT Yea 
444Q@ Xlabel$=#"011] Percent” 
4444 JF Irt\#@ THEN 


4446 Ylabel$S="h/ne% ” 
4457 ESE 

4456 Ylabel€="h/nsmooth” 
446@ END IF 


4464 PRINT “SF 1.5,2:PU PA SO@,-12@ CP” .-LEN(Xlabel?)/2.°@;LE”  Xlabel€:“” 
4466 PRINT “PA -11,5@ CF @,°;,-LEN(Ylebel€)/2eS/6:"D] @,1i LE": VYlebel$S,"* 
4472 PRINT "CP @,@° 

4476 Ipn=@ 

448@ BEEF 

4464 JNPUT “WANT TO FLOT DATA FROM & FILE (1]=5,0@=N)?" Cin 

44P& IcreQ 

4452 IF Otp=1 THEN 

E295 ESeEr 

45@Q INPUT "ENTEF THE NAME OF THE DATA FILE” ,O_files 

4504 BEEP 
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4SCec 
Aci? 
4S1€ 
4522 
4506 
4576 
4S32 
453E 
4542 
4544 
4546 
4SS2 
4956 
4562 
4564 
4568 
4572 
4576 
4582 
4584 
4588 
A292 
4596 
4622 
4604 
4628 
4612 
4616 
4622 
4624 
4628 
4632 
4636 
4642 
4644 
4648 
46S2 
4€S6 
4662 
4664 
4668 
4672 
4676 
4682 
4684 
468¢e 
AES. 
4E9E 
4722 
4704 
4708 
TA aed 
4716 
4720 
4724 
Lact 
ie 
473€ 
4720 
E766 


ONE US  OSELECY (C= INE@F. TecOGwany. , [ios 
ASSIGN @File TC C_ files 
Seer 
INPUT “ENTEP THE BEGINNING RUN NUMBER” ,Md 
SEEF 
INPUT “ENTEF THE NUMBEF OF y¥Y-Y PAIRS STOREO" .Noairse 
Beer 
Piewim ENTER OESTREG HEAT FLU¥Y- 0 
BEEF 
PRINTEF IS 1 
PRINT USING "4x,""Select @ symbol “"' 
PRINT USING “4x,""1 Star Zo Pliue sion” 
PRINTVUSING “4X%,°°3 Circle 4 Squere’’™ 
PRINT USING "4X,""S Rombus”"" 
PRINT USING “4xX,""6& Right-side-up triangle””" 
PRINT USING “4X,""7 Up-side-down triangle””® 
INPUT Sym 
PRINTER 15 785 
PRINT “PU DI” 
IF Sym=1 THEN PRINT “SMe 
TF Sym=2 THEN PRINT °SM+* 
IF Sym=2 THEN PRINT “SMo" 
Nn=wé 
IF Ilog=1 THEN Nn#}] 
IF Md?! THEN 
FOR Jel (0 (Md=1) 
ENTEP @File.Xa,Ya 
NEXT I 
ENO IF 
Q1=Q 
IF llog=1 THEN Q=*L0G6(Q) 
FOR I=] TO Npairs 
ENTER @File:ixa,8(+) 
YaeE(?) 
FOR K*1 TO Nn 
YaetYa+8(kK )#Q*K 
NEXT K 
IF llog=)] THEN Ye=fXPC( Ya) 
IF ITlog*@ THEN Ya=Q1/Ya 
IF Irt#@ THEN 
IF Xa=*@ THEN 
YouYa 
Ya*} 
SUSE 
YaFYa/Yo 
ENO IF 
Scsia 
Hsm=FNHsmocthn(C,Xa,isat ) 
Ya=Ya/Hsm 
END JF 
Xx(I-1 )#Xa 
vl eye 
k=(Xa-XmindeSfa 
Y="VYe-Ymin)«Sfy 
IF Sym-2 THEN PRINT “SM" 
Wee ovine Hen PRINI “SP 1.2.7.4" 
a | Neate oe Yo 
gee WHEN PRINT “SR 1.2.1.6 
Peeovm@e THEN PRINT “UCD .4,.99,0,-6,-4,2.0.8.4,C.. 
eee = Phen PRINT “UC2,¢€,9°.=-23,-6.-2,6.20,6.5."6. 


4 pot 
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Zo 4e 
Cee 
4756 
4762 
4764 
476E 


= oney 
ae 


4776 
4782 
£784 
4788 
47932 
479596 
4822 
4804 
42?@8 
ABZ 
4816 
4€20 
4824 
4828 
4832 
4236 
4642 
4844 
4848 
4852 
4856 
4862 
4864 
4868 
472 
4876 
4882 
4884 
4888 
4e92 
4eg¢ 
£922 
4904 
4906 
4912 
4916 
4920 
4924 
4928 
4°32 
4936 
4942 
4944 
4948 
495. 
4°56 
4960 
£964 
4£96@ 
4972 
4976 
£980 
4954 


IF Sym=£ THEN PRINT “UCC,S.2,99,3,-€,-6.2,2,8. 

IF Sym=7 THEN PRINT “UCC.-S.2,99,-5,65,6,2.-2,-6.° 
NEYT I 

BEEF 

ASSIGN @File TO » 

END IF 

PRINT “PU SM” 

BEEP 


INPUT “WANT TC PLOT @ PDLYNDMIAL (1#1,Q@=N)?" Dic 
IF Gio#i THEN 
Beer 
INPUT “SELECT (Q@=LINEAR,1*#LOG(X,Y))",llog 
Iprnt#l 
CALL Foly(Itn) 
FOR Xee#xXmin TD Xma- STEF xstep/2S 
Ienelentl] 
Ye=FNPoly( Xe) 
Ye  Ya-YmindeSfy 
n=(karXmin)eSf>» 
IF Y<@ THEN Ye? 
IF Y>1@@ THEN GOTD 4868 
Pu=Q 
IF Jpn#=l1 THEN Idf\=Icn MOO 2 
IF Ipn#=2 THEN Idf#Ien MDD 4 
IF Ipn#=3 THEN Idf#Ien MOO 8 
IF Ipn#=4 THEN Idf#Ien MOD 16 
IF Ipn=S THEN Idf*lIen MDD 32 
IF Idf#) THEN Fuel 
IF Pu=@ THEN PRINT °PA*,X,Y, "PD" 
IF Pus] THEN PRINT "“PA",X,Y, "PU" 
NEXT Xa 
PRINT “PU* 
lpn=Ipn+l 
GOTD 4482 
ENO IF 
BEEP 
INPUT ("WANT 710 OUIT )2yY G=N 07> lout 
IF Iquit#i THEN 49@4 


GOTO 4482 
PRINT “PU S®Q" 
SUBEND 

SUB Stats 
PRINTER IS 791 
J=Q 

b=? 

BEEP 


INPUT “PLOT FILE TC ANALYZE?" Files 
ASSIGN @File TO Files 

BEEP 

INPUT “LAST RUN No?/@=QUIT)*,Nn 
IF Nr=@ THEN S292 

NneNn-J 

S$x=2 

Sy=@ 

S7=@ 

Sxs=@ 

Sys=2 

S2 ee 

FDrF Je} TCO Wr 

JeJ+) 
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498E 
499° 
499€ 
See 
5004 
Se0eB 
S@12 
S@16 
5020 
Sec4 
5028 
| lap 
S036 
S24e@ 
S044 
S@48 
S52 
S@s6é 
S@6e@ 
SQ64 
S@6e 
Se@72 
S276 


5282 
SQB4 


ENTEF @File.C,1 

H=Q’T 

SaaS +U 

Sxse57s+0°C 

Sy#Sy+T 

Sye=Sys+T oe 

S2"S 2th 

S2seSostH . 

NEXT ] 

Qave#S./Nn 

Tave#Sy/Nn 

Have®S=-/Nn 

Sdeva=SQP( ABS( (NneS 22-55% °> )/ (Ne (NnW-1]))): 
Sdevt=SOR( ABS( (NneSys-Sy “2 )/(Nne(Nn-] )))) 
Sdevh#=SQR(ABS( (Nn*S28-S272 )/(Nne(Nn-1)))) 
Sh=10@@*Sdevn/Have 

Sq210Q*Sdeva/Qave 

St=1@Q*Saevt/Tave 

IF K=] THEN S@B4 

PRINT 

PRINT USING “1J1X,°°DATA FILE: °° ,14A° i Files 
PRINT 

PRINT USING °“11X,°"RUN Htube SdevH Qdp SdevQ Thetet Sdevl-- 


K= ] 
PRN elnGw) 1x, DO,202),0.50E,1X,30.20),2X,00.30,1%,30.20°: J ,Have,Sh, Gave 


Poo, lave ,St 


52@B8 
Se@s2 
S236 
5182 
5104 
5108 
S112 
S116 
5120 
S124 
Size 
S132 
S136 
5142 
5144 
5148 
5152 
S1S6 
5162 
S164 
Sree 
$172 
5:76 
518e 
S184! 


GOTO 4942 

ASSIGN @Filel TO « 

PRINTER IS } 

SUBEND 

SUE Coef 

Wome Caly7ee.1O.1C).C(10),8(4),N,lornt.000, llog 
Peer 

INPUT “GIVE A NAME FOR CROSS-PLOT FILE” ,Coft 
BEEP 

PIeUleaoolr Ul IYre KO8G vs Ot, l=h vs O02, 2eh vs Q)°,Oac 
CREATE BOAT Cof$,6 

ASSIEN BFile TO Cofs 


BEEP 

Teele SELECT (O=cINEAR, leLOG(X,Y))> ,llog 
BEEP 

TNEUI en eRe OIE PERCENT (-1=STOP)° Bop 
BEEF 


INPUT “ENTER TUBE NUMBER (1, 2, 3, 4, OR S)",JItn 
1F 60p‘°@ THEN S176 

Sell Poiy« lin? 

OUTPUT OF ile. Bop,B/e) 

G6O0TO $144 

ASSIGN @File TG » 

SUBEND 


SIBEMKXXXXXXXXXXXXXXKXXKKMMRXKKX KA KKK KR MK KK XXKXKXKKXKKKKXKKKKX MAK KKK KH 


Siec! 
SI3E 
520? 
5204 
S206 


Eotn 


—_e + & 


S216 


SUB Flot 

Wore tomvee. .e,.1C 70°10) ,€(4),Noc,lecnt Ope, leo 
GiM v.17) 

INTEGEF 11 

PRINTEF IC 1 

BEEP 
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Bll 
ne 
Spe al 
asf 
ee ee 


S200 
5. 4¢ 
6244 
524s 
See 
a Gis 
5260 
5°64 
5.266 


a045 
Poy er 


SL7E 
S2B@ 
5284 
5-88 
6232 
S230 
5302 
5304 
S306 
Bale 
$216 
320 
5324 
S328 
S336 
5242 
$344 
S34E 
e2$2 
S256 
5360 
S364 
$363 
G27. 
S37eE 
S282 
6384 
S28e 
5392 
Se |e 
S402 
S404 
5488 
C412 
S416 
S420 
6424 
S478 
S452 
S426 
5442 
S4éé4 
S4ck 
S452 
S4SE 


Idve] 


INPUT “LIKE DEFAULT VALUES FOF PLOT ¢ti=y1 DEFAULT I.¢=h 


Opoes0 
BEEP 
PRINT USING "4e, 7 7Sélect Gpriem "=~ 
PRINT USING °6%,°°@ ao versus delta-T( DEFAULT )""” 
PPINT USING "6X,°"] hh versus delte-T°"” 
PRINT USING "6X, ° 20 °h vecciicug: — 
INPUT Opo 
EEE 
INPUT “SELECT UNITS (@=S1‘DEFAULT),1"ENGLISH)” ,Jun 
PRINTER IS 7Q5 
IF Idv<21 THEN 
BEEP 
INPUT “ENTEF NUMBEP OF CYCLES FOR X-AXIS",Ca 
Beer 
INPUT “ENTER NUMBER OF CYCLES FOR Y-AXIS" ,Cy 
BEEP 
INPUT "“ENTEP MIN X-VALUE (MULTIPLE OF 10)°,Xmin 
BEEP 
INPUT “ENTER MIN Y-VALUE (MULTIPLE OF 10)°,Ymin 
EESE 
IF Opo=@ THEN 
Cy=2 
Cree 
Xmin=) 
Ymin™]00@ 
END IF 
IF Opo=)] THEN 
Cy#2 
Cx=2 
Kmin=] 
Ymin=10@ 
ENE ae 
IF Opo=2 THEN 
Cy=2 
Cre? 
Xmin=]1000 
Ymin=10Q 
ENO LE 
ENC IF 
BEEF 
PRINT “INS SPIIG IP 2300, 2200, 63500 EbO0..- 
PRINT “SC @.100,@,10@:TL 2;0:- 
Sf>=1Q0Q@/C» 
Sfv=l1@e@/Cy 
BEEP 
INPUT “WANT TO BY-PASS CAGE (1#Y, @=NO( DEFAULT)" ,Ibp 
IF Ibpe!l THEN S9Q@E 
PPINT “PU @,@ PD" 
Nn=9 
FOR J=1 TO Cp+1 
hat=Xminel eo y-1) 
TF JeCy+! THEN Nn=l 
FOF Jz] TG Nn 
IF J=i; THEN PEINT “T 
IF Jer THEM SFINT “°T 
Ae=Xated 
XELOT(Ka/xXmindetf> 
Pett, “PA vo” 0; San 


as 


fa atk 
— —) 
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ar 


5<6C 
S464 
S46€ 
S472 
S476 
S460 
S464 
S486 
Soa 
S495 
SS50¢ 
$504 
5508 
Sain 
216 
$522 
5524 
SS25 
SS32 
S255 
oS4e 
C544 
5548 
SSS2 
SSS6 
ss6e 
SS64 
SS68 
SS72 
SS76 
5580 
S584 
$588 
SE3Z 
6S96 
S622 
S604 
S6QE 
45) @ 
SE16 
$622 
56024 
S62 

S632 
Bo55 
S642 
5644 
S64e 
SAE 
S656 
S660 
$664 
SS6E 
Sere 
Sie Ve 
215219 
So64 
SE8e 
So 3. 
Se oc 


NEXT o 
NEXT ] 
Bruni “PA JOC.C. PU.“ 
PRIN: “PU P& 0,@ PO" 
Nn=S 
Or t#i 1G Cyt] 
Yatz=YminelO°( 1-1) 
IF [=Cv+1 THEN Nn#] 
poem) TO. Nn 
IF J=1 THEN PRINT “TL C 


e* 


IF Je2 THEN PRINT “TL 1 Q" 


YezrateJ 
YeLGT(Ya/Ymin)d*Sfy 
Been orn ey aie tyes 
NEXT J 
NEXT I] 
PRINT "PA @,100 TL @ 2" 
Nn#=9 
FOR I#1 TO Cx] 
XatexminelO*( 1-1 ) 
IF J=Cx+1 THEN Nn#1 
FOR J=#=1 TO Nn 
IF J#l1 THEN PRINT “TL @ 


ce 
< 


tees ei THEN PRINT “TL 0 1° 


Xaexated 
X©L6T(Xa/Xmin)*Sfy 
PRINT “PA°;X,°,100; XT” 
NEXT J 
NEXT ] 


SeINT “PA 100,100 PU PA 100.8 PO” 


Nn=8 
Bor [=] 10 Cy+1 
YateYminel@° (1-1) 
IF IeCy+1 THEN Nn} 
FOR J=1 TO Nn 
IF J=1 THEN PRINT “TL @ 
IF J>! THEN PRINT “TL @ 
YazYateJ 
YeLGT(Ya/YminoeSfy 
PRINT "PD PA 10@,",Y,"Y 
NEXT J 
NEXT |] 
PRINT “PA 100,10@ PU” 
PEIN} “PA @,-2 SR 1.S,2° 
Ti=lLGTt( Xmin) 
ede =1 1G Cx+1 
Xaexminel@*(I-1) 
XeLGT(Xa/Xmin)¢Sfx 
PRINT “PA Ux...” 
IF ITi>=@ THEN PRINT “CP -2, 
IF 1:<@ THEN PRINT “CP -2,- 
Pv=1i+1 
NEXT I 
PPINT “PU FA C,2@° 
Leet ( Y¥min) 
Y1@e1e 
~O" I=} 1G Cyv+i 
YeztYyminel@ Cl-1) 
YeLGT(Ya/YminseSfy 
PRINT “PA 0.°:Y,°" 


2 
| 


le 


7 BAG; PR ae ory) 2 gb OS Te 


yas 


SOLO Vee eee eB loa, 


Cas 
- 


57 0c PRINT "CP =6 >. 2S .xSEr ie Be eee et 

S704 Pel +4 

ST@CE NEYT I 

G7 en CEeEr 

S716) Idlel 

ST2@ INPUT “WANT USE DEFAULT LABELS (1l*eY(DEFAULT),@=N)?”" Id 

Grose JE di HEN 

yee S BEEP 

Bioe INPUT “ENTER X-LABEL” ,Xlabel€ 

S73 BEEP 

5742 INPUT “ENTER Y-LABEL”,Ylabel$ 

744 ENC Ir 

5746 “IF po<c THEN 

§7S2 PRINT “Shot. 2. PUSPA 405-14. 

5756 PRINT “LBC(T: PR -1.6,2 PD PR 1.2,.@ PU: PF .5,-4,.LBwoi PR .S,1; - 

5762 PRINT “LE-T;PR .&,-1:L6sai:PR .5,1.° 

5764 IF Jun#@ THEN 

S768 PRINT “LE) / (K)° 

S772 EESE 

S77 PRINT “LB de/eGR)- 

5782 ENO IF 

$784 ENO IF 

S788 IF Opo=#2 THEN 

6732 IF Iun#=@ THEN 

5796 PRINT "SR 1.5,2;PU PA 40,-46LBq / (W/m:SR 1,1.5;PR €.5,17;3L82;SF 
eo.e3FR O.S7=), Eee 

S822 ELSE 

5804 PRINT °SR 1.8,2:;PU PA 34,-14:LBq / (Btu/hr:PR .5,.5:LE.;PR .§,-.5: 

580e PRINT “LBft;.PR .8.1°SR 1.1.5; LEZ (SR S22. moe see oe 

SEi2 END iF 

S816 ENO IF 

S820 IF Opo#2 THEN 

5824 IF Iun#=@ THEN 

S826 PRINT “SR 1.5,2:PU PA -12,4@8.D1 OQ, 1;:L8q 7 (W/mePR —-190-5. See ome 

B2aSR 1.5.23 Pre. S768)" 

S832 ESE 

SESE PRINT “SR 1.5,2;FPU PA -12,32.01 @.1.UBa 7 CBtUG/he, Pee--=- 5 Cee 
Sa. Su * 

5642 PRINT “LEft; Sk 1.1.5; PR —1. 75262) Pel oor eee 

S644 ENO IF 

S84€ END IF 

5652 IF Opo>@ THEN 

SESe IF Iun=@ THEN 

SS6e2 PRINT "SR 1.5,2;PU PA =<12,3€:01 @. 1 EB 7 Clim. Fea) Ss ee 

SSR Veo re 2 oe. on.” 

Se@64 PRINT “L&..PR .S,@.LEr )” 

Ses8 ELSE 

Sé72 PRINT “SR 1.5,2:PU PA —i2,28:01 @,1.LBh 7 (8tu/he, Pke—-S2 See 
Bt ay Sais 

SE76 PRINT “LBft; PR -1.,.5,5R91, 1.55 LBo. SR Sy ese 

LE)” 

Seee END IF 

Se4 ENO IF 

Sste IF Idl=@ THEN 

Segc PRINT “SF !.£,2-PU Fe SCV=16 CP". =<LEN Sdatbeleiyen G.ce je iene 2. 

SEge PRINT “Pl —14.50 CF 0." rHENCY labels  seS7e. Die eee te 

5S2¢ PRINCE C0" LI" 

£904 ENC IF 

SS9QE Ipn=2 


2 


SSicl Repeat | 


SIIE 
o37e 
S9r4 
S203 
S952 
S936 
S940 
S944 
S948 
Scie 
S956 
SSE@ 
S964 
SSIES 
Slr 
S976 
S982 
S984 
S988 
5392 
S396 
6222 
6204 
6228 
6212 
6@16 
6222 
B@24 
BQ28 
6232 
BQ36 
6@42 
6044 
6248 
Be@S2 
6@SE 
BQE@ 
6Q@64 
B@SE 
6@72 
6277 
6282 
B2@84 
BQ8eé 
B@92 
6@9E 
B1@e@ 
6104 
B1@€ 
Blie 
B1l6 
E122 
6124 
G12E 
Bales 
ESE 
£142 
Gideé 
El4E 


Wiles .ese 
Xul=-!.E+6 
Len2e 

BEEP 

Of =] 


INPUT “WANT TG PLOT DATA FROM & FILE (18Y: OEF AULT), Q@=N)?" ,O} 
IF Obe#!l THEN 

BEEF 

INPUT “ENTER THE NAME OF THE DATA FILE" Of: leSi@) 

ASSIGN @File TO Ofile$(@) 

BEEP 

Npairsecl2 

INPUT “ENTEF THE NUMBEF OF X-Y PAIRS STORED(DEFAULT=2@)”" .Npairs 
BEEP 


Itneltnt+l 
INPUT “ENTER TUBE NUMBER (1, 2, 3, 4, OF S)“,I1tn 
BEEP 


PRINTER IS lI 

INPUT “WANT DEFAULT SYMBOLS? (YES=@ (OEFAULT),NO=] )",Symb 

Sym=lItn+2 

IF Symb#@ THEN 
60TO 6@36 

ENO IF 

PRINT USING “4X,""“Select a symbo]:"“"° 

PRINT USING “&xX,-71 Ster < Plus sign"*” 

PRINT USING “6X,°°3 Circle 4 Squere*”” 

PRINT USING “6X,°°S Rombus““°” 

PRINT USING °6X,"°6 Right-side-up triangle””” 

PRINT USING “6X,°"7 Up-side-down triangle””” 

INPUT Sym 

PRINTER IS 725 

PRINT “PU OI" 

IF Sym=1 THEN PRINT “SMe” 

IF Sym=2 THEN PRINT “SM+" 

LES SVm= So PHEN PRINT. «oto - 

FOF I=] TC Npairs 

ENTEF @File:Xxy(e) 

IF Opo=#@ THEN 
Ya=xy(Itn-]) 
Keexy(ll¢+lItn) 

ENO IF 

IF Qpo#] THEN 
ToeAy Coe) 
Xe=zxy(ll¢lin) 

ENO IF 

IF Opo=2 THEN 
Ya=Xy(S+Itn) 
Xaexy(Itn-1 ) 

END IF 

IF Xac<xll THEN Xll=Xe 

IF Xe>Xul THEN xXul=Xe 

IF ITun=1 THEN 
TF O50.2 THEN Xaedae!. 
IF Qpo+@ THEN YeFrae.} 
IF Ope=C THEN veeYer. 3 
LE Opoe=2 THEN Aa=rvae.3 

ENG sir 

YEeLOT( Xa ’Xmin)eCf> 

YeLor(Ya/YmindeSfy 


73 


E:Sc byee 

G {SomIb ALIA Gumtinmee cr lc Ua 

B16e 60TO B2l2 

E164 IF Sym>Z THEN PRINT “SM" 

BI6B IF Sym<4 THEN FRINT “SR 1.4,2.4° 
G17. U1F lel=O i Her 


B17E Pe UN ee eee 
Seucie Sey a 
61B4 PeClINH@e ee fox igor al 


6168 END IF 
51S8% iF SymeoeTHEN PRIN, SR lee tac 

6196 1F Symed THEN PRINT “UCD, 4.99.00-6--4 0 ce 4 o.° 
6200 IF Sym=S THEN PRINT “UC3,0,99,-2,-6,-2,.6,2,6,3,-6." 
6204 IF Sym=€ THEN PRINT “UC@,S.3,99,2,-8,-6.2,3.6)" 
6208 IF Syme? THEN PRINT “UCQ,-S.3,99,-2,6,6,2,-3,-8:" 
6212 NEXT 1 
Ses oh oie 

220 BEEP 
6224 llabel 
6228 INPUT “WANT TO LABEL? (1=©Y(DEFAULT),Q=N)*,Ilab 

6232 IF llabe\l THEN 


6236 PRINT °SP@.SP2° 
6242 BEEP 
6244 IF Klabe@ THEN 
24B Xlab=ES 
6252 Ylab=€S 
6256 INPUT “ENTER INITIAL X,Y LOCATIONS" ,Xlab,Ylab 
6260 XtteXlabp-S 
6264 Ytte"#Ylab+8 
626e PRINT "SP 1,1.5° 
bere PRINT “SM,PA" ,Xtt,Ytt,°L6é Tube % File” 
6276 Vide yit-5 
6 2B PRINT “PA xtt,Ytt, “LE No O:1l Name” 
62B4 IF Sym]! THEN PRINT “SMe 
628s IF Sym=2 THEN PRINT °SM+”" 
sped Pe IF Sym=2 THEN PRINT °SMo” 
6296 Klab=1 
6302 ENO IF 
E204 K y=] 
E308 CALE Symot(h las.) fab. oym,lelars? 
oA PPINT “SF 1,1.5;5m" 
BZ1E IF Sym<4 THEN PRINT “PR 2,0" 
22 PRINT a FRoe. Gao celta ae 
6224 BEEP 
E32E INPUT “ENTER BOP( @=DEFAULT )” ,Bop 
Bos IF Bop’ 10 THEN PRINT “PR 3,@:LE";:Bop.”” 
£336 1F BSo0p>9 THEN PRINT “PP 1.5,@.LE":Bop:"** 
6242 PRINT “Pee eee Le Oli lecego. * 
6344 PRINT, ~SeO-SPI:SK 155,2- 
ED4E Ylab=Ylapb-S 


E252 Mie NOe 1 Ee 

Sects Sala 

EZ6@ ASSIGN @File 7C ° 
63640 Mill=s pe 

EZ66 Aulseteieese. 

62°< ' SOTC BG460 

Eo7e ENGwie 

E[8C PRINT “PU SM" 
E264 PEEP 

cSBE Gc_on-l 


174 


EZS2. INPUT “WANT TO FLOT A POLYNOMIAL (1=1'° DEFAULT), @EN >? .Gc_on 
E296 IF Gc_on=i THEN 


64090 Beer 

64204 FRINTEF IS ] 

B4QE INPUT “WANT OEFAULT LINE TYFE? « YES*@ (OEFAULT),NG#l>“ Ln 
6412 Ipn=Itn 

B41 IF Lno=@ THEN 

6422 BOTC 6448 

6404 ENG IF 

B47e PRINI@USING “4X,~ "Select line type 
Bao. PRINT USING "6X," "@ Solid line“*" 
6436 PRINT USING “6x,°°1 Dashec”**" 
6440 PERINTTUSING “6x, "<...5 Longer line — dash ““ 
6444 INPUT Ipn 

6448 PRINTER IS 7@5 

6452 BEEP 

6456 Ilog=1 

6462 INPUT “SELECT (@=LIN,1*LOG( DEFAULT) )",Ilog 
6464 Iprnt=l 

6468 CALL Poly(Ofile$(¢),Npairs,Itn) 

6472 FOR xXx#@ TO Cx STEP Cx/220 

6476 XaexXminela°xx 

6482 IF Xa<xlI OR Xa>xXul THEN 6572 
6484 len*Icn+1 

648e PusQ 

6492 IF Ipn=1 THEN Idf#Ien MOOD 2 

6496 IF Ipmel THEN IdfsIcnm MOO 4 

6522 IF Ipm=2Z THEN Idf#Icn MOO 8 

6504 IF Ipmn#=4 THEN Idf#Icn MOO 16 
65@8 IF Ipn*S THEN Idf#Icn MOD 28 
6512 IF Idf=} THEN Pue} 

6516 YasFNPoly( Xa) 

E522 IF Ya<Y¥min THEN 6572 

6524 IF Ium=]} THEN 

Boece IF Opo<2 THEN Xe#=Xae!.8 

6532 IF Opo>@ THEN Ya=Yae.1761 
€S36 IF Opo#@ THEN YaeYae.317 

654@ IF Opo=2 THEN XaexXa*.317 

6544 ENO IF 

6548 YeLBT(Ya/Ymin)*Sfy 

6SS2 X@LBT(Xa/xXmin)*Sfa 

BS5E IF Y<@ THEN Y#@ 

6562 IF Y>1@@ THEN GOTC 6572 

6564 IF Pus@ THEN PRINT “PA*,X,Y, "PD" 
6568 IF Fuel THEN PRINT “PA",X,Y, “PU” 
ES7c NEXT X> 

E576 PRINT “PU” 

6S8@ ENO IF 


poe 4 EESP 

SsscmeNPUT “WANT TC QUIT (1=),Q@=N(OEFAULT ))" Tat 
Soo. oF Iqgt=] THEN 6602 

Bese GOTO 5916 

BeOUsePRINT “PUl PA @,@ SPO” 

6604 SUBEND 

BEQE! 

BES ZIXXXXXXXKYKKXXKXXKM KM MMM MMM KK KIM KKM KKM MMMM MMH KKK MK KK HAIN KD! 
EE16! 

BeoUmeeUe SymbtxX,Y,Svm,Je},t a) 

Beol mr Syn 2 THEN PRINT “SM” 

PeccmmiteoyMce THON PPINT “Sh 1,4,.,4° 


7S 


GEZ2 Yaaee 
6636 1F tj=)] THEN \ac=.8 
6640 1F Icl=@ THEN 


6644 PRIN] “eA oes, to eG, 
65S5 “ELSE 

BBSs PRINT: PA”. &, Ye Yed,. Po” 
Beso sENDeIE 


B66@ IF Sym>3S THEN PRINT “SF 1.2,1.6" 

E664 IF Syme4 THEN PRINT “UCZ.4.96,0,-€,-4.20,C,6,4,0." 
6686 IF SymaS THEN PRINT “UC2 (0 355-5 (-6  —-5 ease ee oe 
6672 IF Syme6 THEN PRIN? “UCO S22, 55.5. —¢ 5 0e jee 
66576 IF Syme? THEN PRINT “UCO{=S 4, 99. -3 276.0, -5,-c 
665Q IF Kj=#]1 THEN PRINT °SM:.PR @,-.8° 

6684 SUBEND 

6688! 

BESDIKXKKXXXKXXXKXKKXK KKK KK KKK KK KKK KKK KKK MM KKK KK KK KKK KKK KKK KKK KM KKK KK KKK KKKK 
6696! 

6700 SUE Fixup 

6704! FILE: FIXUP 


6708 | 
6712 O1IM Emf( 347 7Amp(11),Volt(4) .Leteta) 
6716 “BEEP 


6720 INPUT “BLOSFILE 10 FIXUP (O2f1 les 

6724 ASSIGN ®@File2 TO O2files 

6728 Difiles#® TEST~ 

6722 CREATE BDAT Olfile$,62 

6736 ASSIGN @Filel TD Olfiles 

6740 ENTER ®File2:Nrun,DateS,Ldtc(e) ,Itt,8or,Nht ,Natp ,Nrt ,Corr 
6744 DUTPUT @FileliNrun,DateS,Ldtc(e) ltt, Boo, Nht,Natp,Nrt,Corr 
6748 FOP I=1 TO Nrun 


6752 ENTER @File2.Told&,Emfie) Volt(¢),Amp(e ) 
6756 IF l=) THEN 6764 
6762 OUTPUT @Filel.8op,Told&,Emfie) Volt) Ample) 


6764 NEXT I 

6768 ASSIGN @File2 TO « 

6772 ASSIGN @File! TO ¢ 

6776! RENAME TES!” T6° D2_files 
678@ BEEP 700¢,. 
6784 BEEF 40202,. 
6788 BEEF 4002,. 
6792 SUSEND 
6796! 
BBOO!XXXKXKXKXKXKXKKKKXKKKKKKKK KKK KKK KK KK KK KK KKK KKK KK KKK KK KK KKK KKKKKKKKKKKKKKKK 
68@4! 

6808 SUE Move 

Sel. FP LEe INA rove 

6816 
662@ Olm A(66),8(66:},C(66.,0(66:,E/6E),F(66),6(66),H(66),1¢666),#(66),L°(66 7), m( 66 
) 

6824 DIM N(66) Emf(54), Volt(2) Amal ll d Ldtet4) 

68°68 BEEP 

Gece Wie die CURE LEE TC MOVE spe es 

BESS “GSelGhe eagles 10 Ul _ fs ves. 

S40 ENTER Ori lec; Neer.) olds betes 41tt. bes wnat Nete ieee cae 

6644 FOF lei TC Nerur 


Fy PI Ad 


BS4E ENTEF @Filec.Told$ 

Ges. ENTEP OF alec, All ),6(17,Cl* 01) ECL? File at ee ee 
LNT 

BESE ENTEF @FileciEmfle: Volt’ *i,Ample) 


6€6@ NEXT 1 
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meee o= of. C§ i.e. TO « 


Beoe BEES 
Boge) INPUT “SHIFT O1Sk AND HIT CONTINUE” ,DI 
Soro EEE 


688C INPUT “INPUT BDAT SIZE” ,Size 

6854 CREATE EDAT Do_file$,Sice 

688E ASSIGN @Filel TO Do_files 

See OUTEU) @Filel«Nrunm, Gates Ldtc(*). 1tt,bop .Nnt .Natc,Nrt.Corr 
Secor or t=) TC Neun 


£9200 GUIRUT C§atel: icilds 

6904 OUTEUIew ei leleat iy. 6! 1 ert) 641). FC1),GCl HO). J0C1) KIT), Ll) eM 
Cle) NI ) 

6928 DUTPUT @Filel:Emfl*) ,Volt(*) Amp le) 


Beles NEXT 1 

6916 ASSIGN @Filel TD « 

692@! RENAME "TEST" TO Do_files 
6924 BEEF 200¢,. 
6926 BEEF 4002,. 
6932 BEEP 4002,. 
6936 PRINT “DATA FILE MDVED’ 

6940 SU8END 

6944! 
BS48!KXKXXXKXKXKXKKXXKKXKKXXKKKXKXK KKK KKKK KKK KKK KKK KK KK KKK KK KK KK KKK KKK KKK KKK KK 
booc 

6356 SUB Purg 

6960 BEEP 

6964 INPUT “ENTER FILE NAME TO BE DELETED’ ,Files 

6968 PURGE Files 

6972 6DTD 6962 

6976 SUBEND 

698eQ! 

E9841 XXXXXXKXKXXKXKXKXKKXKKXKK KKK KKK KKK KKK KK KKKK KKK KKK KK MK KK KKK KKK KKK KKK KKK KK 
£988! 

68982 SUE Comb 

6996! FILE NAME. COME 


PPI 2 


700! 
Geese Oi Emf(s4),Volt(2),fmp(11),Ldtcet4) 
70@& BEEF 


eile UNEUT -OLD FILE TG FIXUP” .D2_ files 

7016 ASSIGN @File2 TO DZ_fiies 

mozo. Dl fi PeS="TEST” 

7024 CREATE BDAT Di_file€,30 

7@28 ASSIGN OFilel TO Dl_files 

7@32 ENTER @File2:Nrun,Date?,Ldic(*),1t1,Boo,Nht .Natp,Nrt,Corr 
7056 IF K=@ THEN DUTPUT @Filel:Nrun,Date$,idtc(+¢),1tt,500,Nht ,Nato,Nrit,Corr 
704@ FOR I*1 TO Neun 

7044 ENTEP @F: lect: 5o0p,Told&,Emf(*),Volt(e),Amp(e) 

7046 OUTPUT @Filel:Bop,ToldS ,Emf(*),Volt(*),Ample) 

7052 NEXT I 

7@5E ASSIGN @=:lec TO « 

meee RENAME “TEST” TO Di_files 

7064 BEEP 400C,.2 

7068 S8EEP 

7272 Ora=l 

7076 INPUT “WANT TO ADS ANOTHEF FILE (1=#Y,Q@=N( defau!ti)7",Ovre 
7@5@ IF Ot a=] THEN 

7QE4 Fe} 

7Q@E@ BEEP 

7@92 INPUT "“CIVE NEW FILE NAME” Nfileg 

7@96 ASSIGN OFilel> TO Nfiles 


pays 


71 GC. GOTC, 70352 

71@5 END 1F 

71Q0S ASSIGN @Filer TO * 

foe. SUBENL 

7116! 
PILAUXKXKKKXXXXXX KKK KKKKKXKXK KKK KKK KKK KK KK KK KX KK KKK KKK KX KKKXKKKKK KKK 
Fie Wa a 

712@ SUE Readplot 

7132 DIM Qdp(S),Htube(S), Thetab(5) 

713E€ PRINTER 15 70} 

7140 INPUT “ENTER FILE NAME’ Files 

7144 INPUT “ENTER THE NUMBER OF DATA PAIRS" ,Nrun 
714€ ASSIGN @Filel TO Files 

7152 FOR 1=1 TO Nrun 


cise ENTER ®8Filel:Qdp(*),Htube(*),Thetab(e) 
7160 PRINT Qdp(e) 

7164 PRINT 

7166 PRINT Htube() 

Wife PRINT 

76 PRINT Thetab(¢) 

7182 PRINT 

7184 PRINT 

7188 NEXT 1 


7192 SUBEND 
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